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Abstract 

The strong-coupling character expansion of lattice models is reanalyzed in 
the perspective of its complete algorithmization. The geometric problem of 
identifying, counting, and grouping together all possible contributions is disen- 
tangled from the group-theoretical problem of weighting them properly. The 
first problem is completely solved for all spin models admitting a character- 
like expansion and for arbitrary lattice connectivity. The second problem is 
reduced to the evaluation of a class of invariant group integrals defined on 
simple graphs. Since these integrals only depend on the global symmetry of 
the model, results obtained for principal chiral models can be used without 
modifications in lattice gauge theories. 

By applying the techniques and results obtained we study the two- 
dimensional principal chiral models on the square and honeycomb lattice. 
These models are a prototype field theory sharing with QCD many prop- 
erties. Strong-coupling expansions for Green's functions are derived up to 
15th and 20th order respectively. Large-A^ and N = oo results are presented 
explicitly. Related papers are devoted to a discussion of the results. 
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I. INTRODUCTION 



It is certainly appropriate to consider two-dimensional principal chiral models as a the- 
oretical physics laboratory. These models display a rich physical structure, and share with 
four-dimensional gauge theories a number of fundamental properties: nonabelian symmetry 
with fields in the matrix representation, asymptotic freedom, dynamical mass generation. 
Moreover, principal chiral models admit a. 1/N expansion and a large- limit which is a 
sum over planar diagrams, in total analogy with nonabelian gauge theories. 

However, the absence of local gauge invariance and the reduced number of dimensions 
make chiral models much simpler to handle both by analytical and by numerical methods. 
Moreover, the on-shell solution of the models is known by Bethe-Ansatz methods: a factor- 
ized S'-matrix exists and the particle spectrum is explicitly known. We can therefore try to 
make progress, both in analytical and in numerical techniques, by testing these methods on 
chiral models and, in case of success, applying them to four-dimensional gauge theories. 

The spirit of this approach is well expressed in the papers by Green and Samuel 
who advocated a systematic study of lattice chiral models as a preliminary step towards an 
understanding of lattice gauge theories, especially in the large- limit. One of the techniques 
favored by the above-mentioned authors was the strong-coupling character expansion. How- 
ever, the existence of a large- iV phase transition from strong to weak coupling phase seemed 
to indicate at that time an obstruction to further pushing this method of investigation. 

In much more recent times, a few facts came to suggest that this "no-go" result might 
be over-pessimistic. It was indeed observed by the present authors that scaling of 

physical observables is present in finite- chiral models already in a coupling region within 
the convergence radius of the strong-coupling expansion. Moreover, a change of variables 
corresponding to adopting the so-called "energy scheme" for the definition of the temperature 
smoothens the lattice /3-function to the point that asymptotic scaling is observed within the 
strong-coupling region. These patterns are unaffected by growing A^, and therefore survive 
the large- A^ phase transition. These "experimental" observations led us to reconsidering the 
possibility that a strong-coupling approach could be turned into a predictive method for the 
evaluation of physical quantities in the neighborhood of the continuum fixed point of the 
models. 

A second theoretical motivation for a renewed effort towards extending strong-coupling 
series of chiral models, especially for large A^, comes in connection with the possibility that 
the above-mentioned transition, while uninteresting for the standard continuum physics, 
may be related to a description of quantum gravity by the so-called "double scaling limit" 
|l7HlO|. In simple models, this limit is studied by analytical techniques, but more complex 
situations might need perturbative methods, and strong-coupling seems well suited for such 
an analysis, which corresponds to exploring the region in the vicinity of the first singularity 
in the complex coupling constant plane. 

Another significant change, of a completely different nature, has occurred since the orig- 
inal studies on the strong-coupling character expansion (cf. Ref. |]TT| for a review) were 



performed. The increased availability of symbolic-manipulation computer programs and the 
enormous increase in performance of computers have now made the strong-coupling expan- 
sion a plausible candidate for an algorithmic implementation, that might extend series well 
beyond the level that can be reached by purely human resources, while granting a definitely 
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higher rehabihty of results. 

The purpose of the present work is to set the stage and make a significant effort towards 
a complete algorithmization of the strong-coupling character expansion. 

Two major classes of problems must be handled and solved. The first has to do with 
counting the multiplicities of terms appearing in the expansion. It is basically a geometrical 
problem and it leads to the definition of a "geometrical factor". We must stress the fact 
that this geometrical factor depends only on the lattice connectivity, and therefore applies 
without any modification to the strong-coupling expansion of all spin models admitting a 
character-like expansion, including 0{N) and CP^~^ models with nearest-neighbour inter- 
actions |]T2[. We have completely solved this problem, with no conceptual restrictions on 
the dimensionality and connectivity of the lattice. We have not addressed the corresponding 
problem for lattice gauge theories, but we are confident that no major conceptual obstruction 
should arise in pursuing that program. 

The second class of problems is related to the evaluation of group integrals that appear 
as coefficients of the expansion. Evaluating group integrals is an algebraic problem, and in 
principle a solved one. However, algorithmic implementation is not in practice a trivial task, 
and therefore we limited ourselves to a general classification and to an explicit evaluation of 
the cases of direct interest to our calculations, with a few useful generalizations. We stress 
that the evaluation of "group-theoretical factors" is universal, and results may be applied 
as they stand to lattice gauge theories. 

The representation of the strong-coupling expansion in terms of explicitly evaluated 
geometrical factors and symbolically denoted group-theoretical factors can be achieved by a 
fully computerized approach, and applies as it stands to all nonlinear sigma models defined 
on group manifolds. This is probably the main result of the present paper. However we shall 
not exhibit here the explicit general formulae resulting from our approach, because they are 
so long that their pratical use does necessarily involve computer manipulation; therefore we 
shall make available our results in form of computer files, publicly available by anonymous 
ftp on the host ftp.difi.unipi.it, in the directory pub/campo/StrongCoupling. 

The application of our results to O(A^) and CP^~^ models is definitely simpler then 
the case discussed here, since the evaluation of group-theoretical factors lends itself to a 
completely algorithmic implementation. The corresponding results will be presented in a 
forthcoming paper. 

The present paper is organized as follows. 

In Sect. |I| we review the character expansion, fix our notation and present some useful 
formulae. 

In Sect. |TT|we outline the procedure of the expansion by identifying the logical steps and 
defining the relevant geometrical and algebraic objects entering the computation. Among 
these we introduce the basic notion of a skeleton diagram, whose multiplicity is the geomet- 
rical factor and whose connected value, or potential, is the group-theoretical factor. 

In Sect. 1^ we explain how one may algorithmically evaluate the geometrical factor. 

In Sect. |V|we introduce the problem of computing the group-theoretical factor. 

Sect. is devoted to some technical remarks on group integration. 

Sect. |V11| offers some details on the computation of potentials for principal chiral models. 

In Sect. [VIII| we analyze the main features of the strong-coupling expansion of the two- 
point fundamental Green's functions, introducing a parametrization for the propagator in 
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the case of a two-dimensional square lattice. 

In Sect. IX we discuss the relevant features of the honeycomb lattice, and we present a 
few results for physical quantities. 

Appendix ^ is devoted to a presentation of some of our results concerning the explicit 
evaluation of potentials. 

Appendix ^ is a list of potentials ordered according to their appearance in the strong- 
coupling expansion. 

Appendix is a presentation our results for large but finite N in the square-lattice 
formulation of the models. 

Appendix clarifies some non-standard features of honeycomb-lattice models using the 
Gaussian model as a guide. 

Appendix ^ is the same as App. ^ for the honeycomb-lattice formulation. 

The present paper is the first of a series of papers devoted to the strong-coupling analysis 
of two-dimensional lattice chiral models. In a second paper we will present our analysis of 
the large-iV strong-coupling series by series-resummation techniques, while a third paper 
will be devoted to a comparison with Monte Carlo studies of the large- critical behavior. 



II. THE CHARACTER EXPANSION: GENERALITIES 

The strong-coupling expansion of field theories involving matrix-valued fields and en- 
joying G X G group symmetry is best performed applying the character expansion, which 
reduces the number of contributions to a given order in the expansion and decouples the 
geometrical counting of configurations from the group-theoretical factor. 



The whole subject is reviewed in detail in Ref. [^, and we recall here only those proper- 
ties that are essential in order to make our presentation as far as possible self-contained. We 
shall only discuss the symmetry groups G = U(A^): extensions to SU(A^) can be achieved 
following Ref. fl]] and applying the results presented in Ref. p. 

In the theory described by the lattice action 

Sl = -Ar/?^[Tr{f/(x) U^ix+fi)} + Tr{t/(x+/i) U^{x)}], (1) 

X,fJ, 

the character expansion is achieved by replacing the Boltzmann factors with their Fourier 
decomposition 

exp{N(3[TT{U{x) U^{x+fi)} + Tr{f/(a;+/i) U^{x)}]} 

= exp{N'F{f3)}J2d^r)ZirmXir)mx)U\x+fi)), (2) 

(r) 

where 

= ^ In/ dUexp{N(3{TTU + TTU^)}= ^ Indet J,_,(2Ar/3) (3) 

is the free energy of the single-matrix model, J2{r) denotes the sum over all irreducible 
representations of G, X{r){U) and d(^r) are the corresponding characters and dimensions 
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respectively, and Ij-i = 1,...,N) is a NxN matrix of modified Bessel functions. We 
recall here the orthogonality relations for representations: 

/ dU!^il\U)^^i*{U) = ^6^rUs)Sa,cS,,a, Xir)iU) = ^M(f/). (4) 

In U(A^) groups (r) is characterized by two sets of decreasing positive integers {/} = 
/i, Is, {m} = mi, rrit and we define the order n of (r) by 

s t 

77, = 714. + 714. = ^ /j, ?7._ = ^ rrij . (5) 

i=i j=i 

We may define the ordered set of integers {A} = Ai, Aat by the relationships 

Afc = /fc, k < s; 

Afc = 0, s<k<N-t + l; 

\k= -niN-k+i, k>N-t+l. (6) 
It is then possible to write down explicit representations of all characters and dimensions: 

^ det|exp{i0^(Aj + A^- , 
^ det|exp{^0,(iV-j)}| ' ^ ^ 

= n.<,(.-.) = «^'<'>' 

where (pi are the eigenvalues of the matrix U. As a consequence, it is possible to evaluate 
explicitly the character coefficients z^r) by the orthogonality relations 



^{A}5{A} = / dUexp{N(3{TTU + TTU^)}xix}{U) exp{-N^F{f3)} 
detIj_i{2N(3) ■ 



(9) 



Eq. (^) becomes rapidly useless with growing A^, due to the difficulty of evaluating 
determinants of large matrices. It is however possible to obtain considerable simplifications, 
in the strong coupling regime and for sufficiently large A^, when we consider representations 
such that n < N. In this case, character coefficients are simply expressed by |l|] 

duMkiM = —^^m^{^}iNPr[l + 0{P''')], (10) 

where we have introduced the quantity cr^i}, the dimension of the representation li, ...,ls of 
the permutation group, enjoying the property 

J X{/}(f/)(Tr f/t)W = . (11) 

It is important to notice that the strong-coupling expansion and the large-A^ limit do 
not commute: large- A^ character coefficients have jumps and singularities at /3 = | and 
therefore the relevant region for a strong-coupling character expansion is /3 < |. 
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A consequence of Eq. (|T0|) is the relationship 

;.^5i;o(/9)=/5 + 0(/32^+i), (12) 
and in turn, because of the property 

^(/?) = ^^W), (13) 

one may obtain the large- iV relationship 

F(/5) =/?2 + o(/?2^+2). (14) 

According to the above observations, at = oo the relationship F = (3'^ may only hold 
when /5 < |, even if the function is perfectly regular for all /?. 

For the purpose of actual computations it is convenient to have expressions in closed 
form for the quantities cT{/} and d^i-m} not involving infinite sums or products even in the 
N ^ oo limit. We found such expressions in the form 

and a similar relationship for a^m}- Notice that these quantities are independent of A^. Now 
by manipulating appropriately Eq. (H) we can show that 

n+! n_! 

where 

_^ {N-t-i + h)\ ^ {N-s-j+m,)\ N + l-z-j + h + m, 

^^'-> - 1\ ^iv^TT-yr l[ (n-s-jV. l[ l[ AT + i-.-, • ^''^ 

The essential feature of Eq. (|1^) is the possibility of extracting results with a finite number 
of operations even in the large- A^ limit. As a byproduct we obtain the large- A^ character 
coefficients in the useful form 

(18) 



III. OUTLINE OF THE PROCEDURE 

The general purpose of the strong-coupling expansion is an evaluation of the Green's 
functions of the model as power series in f3. If we are interested in the mass spectrum of the 
model, we may focus on the class of two-point Green's functions defined by 

G^r){x) = -^{xir){U\x)Um), (19) 
«(r) ^ ' 
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and even more specifically we may decide to evaluate the fundamental two-point Green's 
function 



G(x) = -(Tr([/t(x)f/(0))). 



(20) 



Evaluating such expectation values by the character expansion involves performing the 
group integrals that are generated from choosing an arbitrary representation for each link of 
the lattice. As a consequence of Eq. (p!OD, only a finite number of nontrivial representations 
contribute to any definite order in the series expansion of G{x) in powers oi f3; we must 
however find a systematic way of identifying the relevant contributions. 

As a preliminary condition for the definition of an algorithmic approach to the strong- 
coupling expansion of G{x), it is convenient to identify explicitly all the logical steps of such 
a computation and define a number of objects that play a special role in it. 



Each lattice integration variable U{y) can only appear in the integrand either through 
the representation characters defined on links terminating on the lattice site y or through 
the observable whose expectation value is to be evaluated. According to the rules of group 
integration, nontrivial contributions are obtained only if the product of all representations 
involving U{y) contains the identity (the trivial representation). 

We define an assignment {r} to be a choice of a representation for each link of the lattice 
that is consistent with the above requirement. Necessary conditions for an assignment can 
be obtained by a close examination of the rules for the composition of two irreducible rep- 
resentations of U(iV). When we consider Green's functions in the class defined by Eq. (|19[), 
we recognize that the operator whose expectation value we are evaluating, when considered 
from the point of view of group integration, plays the role of a unit length link connecting 
the sites x and 0, weighted with a factor ci^^. Therefore all the relevant group integrals 
can be put into correspondence with integrals appearing in the character expansion of the 
partition function (possibly in higher dimensions). 

Changing the convention for the orientation of links changes each representation r into 
its conjugate f (/<->• m), but, since Z(^r) = ^(f), it does not affect the expansion. Hence we 
can consider all links terminating in a given site as "ingoing". It is now possible to prove 
that an assignment must satisfy the following conditions at each lattice site: 



where the summation is extended to all ingoing links. 

Order by order in the strong-coupling expansion, the relevant assignments involve non- 
trivial representations only on a finite number of links, which allows the possibility of drawing 
on the lattice the diagram of each assignment. Such a diagram is characterized by vertices, 
where more then two nontrivial representations meet, and paths, i.e. chains of links con- 
necting vertices. Orthogonality of representations implies that the choice of representation 



A. Assignments 




(21) 



n± < ^ (non-backtracking condition) , 



(22) 
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FIG. 1. All the non-trivial disconnections of a diagram. 



along a given path cannot change. We will denote the length of each path p by Lp, and 
the corresponding (nontrivial) representation by r^. The topology of a diagram may be 
represented by the connectivity matrix between its vertices. As we shall show later, the 
value of the group integral associated with each assignment can only be a function of r^, Lp, 
and y-, we shall denote it by Ru-n- 



B. Configurations 

The set (n+, «_) does not in general identify uniquely a representation. It is convenient to 
define oriented configurations: they are the sets of all assignments having the same (n^,_, nL) 
for each link of the lattice. The relevance of oriented configurations in the context of the 
strong-coupling expansion stays in the fact that they are in a one-to-one correspondence 
with the monomials one would obtain in the integrand after the series expansion of the 
Boltzmann factor in powers of (3. They are therefore the simplest objects admitting a 
meaningful definition for their connected contributions. 

Eq. ([T0| ) tells us that the lowest-order contribution of any character coefficient to the 
strong-coupling series depends only on tt, = n+ -|-?7,_. Hence it is useful to define (unoriented) 
configurations by summing up all the oriented configurations characterized by the same value 
of ni for each link of the lattice. The set {n} = (rii, ...) uniquely identifies a configuration. 
We might have defined configurations directly as the sets of all assignments sharing the same 
{n}; our procedure insures us about the possibility of defining the connected contribution 
of a configuration. 

We may introduce the diagrammatic representation of oriented configurations, by draw- 
ing each oriented link {n^, n_) as a bundle of n links, of which n_|_ bear a positively-oriented 
arrow and n_ bear a negatively-oriented arrow. Removing the arrows leaves us with a di- 
agrammatic representation of (unoriented) configurations. One may easily get convinced 
that the algebraic notion of disconnection turns out to coincide with the geometrical one. 
In this representation, a disconnection is a set of subdiagrams such that their superposition 
reproduces the original diagram. An example of disconnection is drawn in Fig. 

Without belaboring on this topic, which is widely discussed in the literature [|13|,|14|, we 
only remind that the connected part of a collection of n (abstract) objects is recursively 
defined by the condition that the set of n objects coincides with the sum of the connected 
parts of all its partitions, including the collection itself. In presence of multiple copies of 
the same object, in standard perturbation theory a combinatorial factor appears, which is 
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hidden in the character expansion; as a consequence, when subtracting disconnections one 
must take care of dividing by the corresponding symmetry factors in order to restore the 
correct normahzation. 

The definitions imply that the geometric features 5^ and {L} of all assignments belonging 
to a given configuration are the same; therefore the path p of a configuration is characterized 
by Lp and by the value of the order of r^. 



C. Skeleton diagrams 

It is convenient to reduce each configuration to its skeleton diagram, whose links are 
the paths of the configuration. The topology ^ is obviously unchanged, and each link is 
characterized by the pair of numbers (n^, Lp). 

In order to clarify the relevance of such a definition, let us consider the problem of 
evaluating the group integrals R^^ for the assignments belonging to a given configuration. 

An elementary consequence of Eq. is the evaluation of the simplest nontrivial group 
integral entering our calculations: 

/ dUxir){A^U) Xir){U^B) = XirM^B). (23) 
By applying repeatedly Eq. along the paths we easily obtain 



<2, = |ni%j'"|4?'. (24) 

where v the number of paths of the configuration (with n 7^ 0), and S^^^^"^ is the value of 
the group integral associated with the skeleton diagram, in which all links are assigned unit 
length and weight, and representations are chosen according to the assignment. Further 
simplification is obtained by noticing that the effective strong-coupling variable in the char- 
acter expansion is z{[3) (for large N actually z{(3) ~ [3 because of Eq. (|12|)). Therefore by 
replacing the character coefficients with the ratios 

= 5^ (25) 

we may express the strong coupling series as a series in powers of z, with coefficients that 
are functions of Z{r)] by the way, these quantities for large enough are pure numbers, 
dependent on N but independent of /3, because of Eq. ([18|) . We can rewrite Eq. (^^ as 



flif2, = --^'-'' n4:, (26) 



since Zi-^ = 1, there is no dependence on the lengths of the links with n = 1, apart from the 
overall factor of 2, depending only on the total length of the configuration L = J2p ^pLp. 
Therefore the corresponding Lp indices can be dropped, thus defining a reduced skeleton. The 
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contribution of a configuration to the functional integral is simply the sum of the contribu- 
tions of all its assignments. It then follows from Eq. (^) that whenever two configurations 
can be reduced to the same skeleton, they will give the same contribution. 

An exchange in the ordering of the vertices will not change the topology of a skeleton 
diagram; therefore configurations that are related by this symmetry will give the same contri- 
bution. Moreover, configurations sharing the same reduced skeleton will give contributions 
differing only by an overall proportionality factor, depending on the total length L. We 
can group together all configurations with the same reduced skeleton (taking into account 
the abovementioned symmetry) and the same value of L: their number is what we call the 
geometrical factor. The common value of each of these configurations is proportional to the 
group-theoretical factor of the reduced skeleton: 

«,= E nk.l'-Si? (27) 

{r} P=l 
{n}fixed 

with a proportionality factor . 

The strong-coupling character expansion of a group integral can therefore be organized 
as a series in the powers of z = Zi;o, with coefficients obtained by taking sums of products 
of geometrical and group-theoretical factors. In order to understand the computational 
simplifications achieved at this stage, let us only notice that, at different orders in the 
expansion, the same reduced skeletons may appear again and again in association with 
different values of L; their group-theoretical factors however are computed once and for 
all, while extracting the geometrical factors is a task that, as we shall show later, can be 
completely automatized. 



D. Superskeletons 

Both for the purpose of bookkeeping and in view of the problem of actually computing 
the group-theoretical factors, at this stage we need a classification and labeling of (reduced) 
skeleton diagrams, which must keep track of their topological properties and try to put 
into evidence whatever further simplification we may conceive. We found it convenient to 
isolate for each topology 5^ a "core" topology ^ which we call superskeleton, defined by 
the condition that each vertex in it is connected by at most a single link to any other vertex 
(i.e. the entries in the connectivity matrix are either or 1). 

The essential ingredient for the reduction of a skeleton to a superskeleton is the extraction 
of bubbles, defined as sets of two links in a skeleton connecting the same pair of vertices. Let 
us now recall the decomposition rule for a product of characters: 

Xir)iU)xis)iU) = j:C^rst)Xit)iU), (28) 
t 

where C(^rst) is a set of integer numbers counting the multiplicity of (t) in the product of 
representations (r) ® (s). For all assignments of (r), (s) consistent with a given skeleton, (t) 
must be such that the triplet (r), (s), (i) satisfies Eqs. (pi]) and (^). Therefore replacing 
a bubble with a single link and allowing for all X(t) obtained from Eq. ( PBD to be inserted 
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FIG. 2. Superskeleton topologies. 



in it defines new consistent assignments. Notice however that in general we may not expect 
all these assignments to belong to the same skeleton, since n may vary within the class of 
admissible {t). 

We can repeat the procedure, replacing paths with links when needed, consistently with 
orthogonality of representations and Eq. (p3|), until all the bubbles in the skeleton have 
disappeared. The resulting diagram is the superskeleton of our original diagram. We must 
stress that a superskeleton is not a skeleton diagram, because it does not make sense to 
assign a value of (n, L) to its links. It is however important to observe that the value 
of the group integral corresponding to any assignment {r} on the skeleton 5^ can be 

expressed as a weighted sum of factors S)^^^^ corresponding to the consistent assignments of 
the superskeleton S/' ^ with weights that are related to the factors 

q..) %^ (29) 

obtained by replacing bubbles with single links. 

A superskeleton is completely identified by its topology, and it is worth mentioning that, 
as in the case of skeletons, supersksletons differing only by a permutation of vertices are 
equivalent, and therefore they can be reduced to a standard form. The number of different 
superskeletons that are relevant to a given order of the strong-coupling expansion is bound 
to grow with the order; however for sufficiently low orders their number is so small that we 
found it convenient to label superskeletons by capital letters, in many cases related to their 
actual shapes. A provisional list of labelings is provided by Fig. ^. 

This is the starting point of our classification scheme for skeletons. Reduced skeletons 
are named by the symbol denoting the topology of their superskeleton; the full information 
concerning superskeleton links, denoted by cr, will appear as arguments; using the pair of 
integers ij to denote the link connecting node i to node j, with node numbering fixed by 
Fig. H, the skeletons will be named 

W(cr), Y(cri2; o^z, cxsi; ctm; 0-24; £^34), 

; (^23; cr34; cr4i; cris; o"25; ctss; <^\^i H((Ti2; (723; CT34; (J41; a\^\ G2h\ cae; <^a&\ cse), 

L(cri2; (723; (T34; (T41; (T36; Cr46; (T15; (T25; (J^q\ 0"45)5 R'(o"12; 021,\ CTsi; (Ju; (T24; C25; (^SS; Cr34; C45)- 
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( p» n,p + g, [ni,Li;n2,L2] •—{ ^ n,p, [rii, Li; [ria, L2; rig, L 

n2,L2 



ni,Li 




FIG. 3. Examples of bubbles. 



a contains information about n; for n 7^ 1, also about the length L in the original skeleton 
and the bubble content. For reasons to be clarified later, we need not consider bubbles along 
n = 1 lines. 

In general, a bubble will be denoted by [o"i; (72], where (Xi and cr2 contains the information 
about the bubble links. In summary, a link information will take one of the forms 

a=l {n = l), (30) 

a = n,L (n 7^ 1, no bubble insertions), (31) 

a = n,L, [(TiXi ^^"1,2] [c"2,i; o"2,2]--- (one or more bubble insertions on a line), (32) 

c" = cri^2] [cr2,i; o"2,2]--- (one or more bubble between two nodes), (33) 

the (Tjj themselves taking one of the above forms; insertion of b identical bubbles will be 
denoted by exponential notation, i.e. [cri;cr2]^. Examples of this notation are illustrated in 
Fig. § 

E. Potentials 

As we mentioned before, the possibility of defining the skeletons as sets of oriented 
configurations insures us about the fact that we may consistently define the connected 
contribution of each skeleton diagram to the vacuum expectation value of an observable. 

Since the geometrical notion of a disconnection only depends on the topology of a dia- 
gram, as a consequence of definitions we can define the (algebraic) connected contribution 
of a skeleton starting from its geometrical formulation. As a matter of fact, it is most 
convenient to exploit the fact that n = 1 lines cannot be split, and define the connected 
contribution of a reduced skeleton, i.e. the connected group-theoretical factor, which we shall 
call potential: 



^{n,L} 



{r} P 
{njfixed 



(34) 



connected 



An example of the chain leading from an assignment to the superskeleton and to the potential 
is illustrated in Fig. ^. 

When we are evaluating the skeletons contributing to the partition function, the sum 
of their potentials with the same geometrical factors is just the free energy. Unfortunately, 
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Unlabeled links are (1;0) 
Assignment 



Oriented configuration 



Configuration 




1,1 




(i;o) 



(i;o) 



Skeleton Superskeleton 

FIG. 4. Steps showing that a sample assignment contributes to the potential 
Y(l, [1; 2, 1]; 1; 2, 4, 1; 1; 1; [1; 2, 2]) = W(2, 1) Y(l; 1; 2, 4, 1; 1; 1; [1; 2, 2]). 



there will be in general no correspondence between the connected contributions to an arbi- 
trary Green's function and the corresponding contributions to the free energy. A notable 
exception is that of the fundamental two-point function G{x). In this case no disconnection 
of the vacuum diagram can split the n = 1 line associated with the fundamental character 
Ti U'^ (x) U{0), and there is therefore a one-to-one correspondence between the connected 
contribution of a given skeleton diagram and the contribution of the associated vacuum 
diagram to the free energy. Moreover the weight di.^ = is the correct normahzation, 

insuring that in the large- limit finite contributions to the Green's functions correspond to 
finite contributions to the free energy. From now on we may therefore focus on the evaluation 
of potentials related to vacuum skeleton diagrams. 

It is worth mentioning that we might have defined oriented potentials, but this notion, 
while conceptually useful, does not find any use in our actual computations. 

A final observation concern notations: we shall label potentials with the same symbols 
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adopted in the labeling of the corresponding reduced skeletons. 

We must draw some attention to the fact that our definition of potentials, although 
referred to unoriented diagrams, is originated by the problem of evaluating Green's functions. 
Therefore we are assuming that the orientation of one of the links has been fixed. By a trivial 
symmetry of conjugate representations, our potentials will be one half of the corresponding 
vacuum contributions to the free energy. 

Including this factor of 2, the disconnections drawn in Fig. ^ can be written as 

disc(2W(2, Li, [2, L2; 2, 0, 1])) = 2 x 2^W{1) W(2, Li, 1) + 2 x 2^W{1) W(2, L2, 1) 
+ 2^W{1) W(2, L1+L2) + I X 2^W{l f. (35) 



IV. COMPUTING THE GEOMETRICAL FACTOR 

The enumeration of all configurations possessing the same reduced skeleton can be com- 
pletely automatized by the following considerations and procedures. 



Eqs. (^) and (^) insure us about the existence of a (non necessarily unique) non- 
backtracking random walk of length J2p ^pLp reproducing the diagrammatic representation 
of each configuration. We therefore generate all non-backtracking random walks with fixed 
length, fixed origin 0, and fixed end x, and we compute the corresponding configuration {n}, 
i.e. we compute n; (the number of times each link is visited) for each link of the lattice. We 
now compare the generated configurations, and discard multiple copies, choosing one (and 
only one) walk for each different configuration. 

The total (bulk) free energy can be computed by summing over all the configurations. 
Therefore the free energy per site can be computed by summing over all the configura- 
tion that are not related by a translation. These are easily obtained by generating all 
non-backtracking closed random walks touching a given site, identifying the corresponding 
configurations, and chosing one configuration for each equivalence class under translation 
symmetry. From this point on, the computation is identical both for the Green's function 
and for the free energy. 

We must notice that at this point we have generated all the sets {n} obeying Eqs. (^I]) 
and (p2[), but not all of them lead to nonvanishing group integrals; we get rid of these "null" 
configurations by defining their group-theoretical factor to be zero. Our computer program 
recognizes and automatically discards two classes of null diagrams: 

Diagrams that can be disconnected by removing a single node. A very simple property 
of invariant group integration allows for the possibility of setting a single integration vari- 
able to 1. As a consequence, one may prove that, whenever the removal of a vertex in a 
skeleton leaves us with disconnected subdiagrams, the value of the group integral factorizes 
into a product of terms that are just the values of its disconnected parts. Therefore, the 
corresponding potentials vanish identically. 

Diagrams that can be disconnected by removing two links, unless the links share the 
same value of n. Such diagrams vanish as a trivial consequence of the orthogonality of 
represent ations . 

Examples of this phenomena are drawn in Fig. ^. 
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FIG. 5. Two null configurations: the first can be disconnected by removing a single node; the 
second can be disconnected by removing a n = 1 and a n = 3 link. 

We compute the reduced skeleton of each of these configurations. We now group together 
all the configurations originating equivalent reduced skeletons (i.e. which are equal apart 
from a permutation of vertices); the geometrical factor is the number of such configurations, 
and we choose one representative configuration for each group. 

We factorize each skeleton "cutting" along n = 1 paths, identify bubbles according to 
the scheme of Subs. |111 D[ and compute the connectivity of the corresponding superskeleton. 
The superskeleton is then either identified as in Fig. |^, or shown to originate from a null 
configuration. Finally we put together all this information to obtain the potential, and use 
the superskeleton symmetry to bring it in a standard form. 

While the data needed in the intermediate stages of this computation can be extremely 
large, the results of the last step (potentials and geometrical factors) are rather compact 
and can be stored for further processing. 

The limiting factor in this procedure is the available RAM. On a workstation with 140 
Mbytes of RAM we were able to generate Green's functions up to 18th order and the free 
energy up to 20th order on the square lattice (they of course involve several new superskele- 
tons beside those drawn in Fig. H). Computer time is not a limiting factor, since the longest 
computations take about one CPU hour on a HP- 730/125. 

At this stage we must clarify what we mean by standard form of a superskeleton. In 
sufficiently complex cases, an ambiguity may arise as a consequence of different sequences of 
elimination of the bubbles. While the resulting superskeleton is always the same, equivalent 
skeletons may receive superficially inequivalent labelings. We have not made an effort to 
reduce completely all these different namings to a standard form, but we were satisfied with 
the reduction to a common form in most cases. We checked explicitly that the computed 
values of differently labeled equivalent potentials are equal. 

V. COMPUTING THE GROUP-THEORETICAL FACTOR 

In contrast with the previous Section, we must say that the evaluation of potentials is 
not yet fully automatized. 

We can routinely generate all the sets {/; m} needed to identify the representations of 
U(A^) to a definite order. The closed formulae presented in Sect. |T| enable us to evaluate 
automatically of their dimensions and their large- A^ character coefficients. 

We can perform the decomposition of the products of these representations, thus identi- 
fying the coefficients C(^rst) and the factors defined in Eq. (pOj). We can therefore reduce the 
evaluation of the group-theoretical factors, by computer manipulations, to a linear combi- 
nation with known coefficients of factors 5*1,^.1, that are nothing but group integrals corre- 
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sponding to consistent assignments of representations on the (unit length, unit weight) hnks 
of a superskeleton with topolc^y S/' . 

Computing the factors is basically a sophisticated exercise in group integration, 

and it is therefore completely solved from a conceptual point of view. The group integration 
over a multiple product of representations can always be performed by decomposing the 
product into sums of representations, via the introduction of appropriate Clebsch-Gordan 
coefficients, and applying orthogonality of representations (Eq. (^) in the last step. This 
may however become a very inconvenient procedure, essentially because of the fantastic 
proliferation of indices (all to be finally contracted, but appearing at intermediate stages 
already in the simplest examples) resulting from writing higher-order representations in the 
basis of polynomials of the fundamental representation. 

We have not seriously tried to overcome this problem in general, i.e. we have no algorithm 
capable of generating the Clebsch-Gordan coefficients for the decomposition of the product 
of two arbitrary representations of U(iV), which would allow to implement the relevant group 
integrations in a computer program. Instead we followed a slightly different approach, more 
limited in purpose and simpler to implement, within our self-imposed limits, without fully 
computerizing the computation. 

The essentials of our approach are the following. 

We observed that, for not very high orders of strong coupling, only a small number of 
superskeletons and low-order representations enter the calculation. Therefore, by making 
use of a few well known results of group integration (that can basically be reported to the 
knowledge of the six- matrix deWit-'t Hooft integral |l^), we managed to compute explicitly 
all the factors S)^^^^ entering in our calculations. 

However, the possibility of inserting bubbles and varying the lengths allows the genera- 
tion of a huge number of different skeletons even starting from a very small set of assignments 
on a superskeleton. The group-theoretical factors of these skeletons can thus be evaluated 
symbolically on wide classes, as functions of the above parameters (which are the same en- 
tering the labeling of skeletons), and the explicit evaluation of the potentials entering an 
actual calculation can be implemented in a computer algebra program. 

The procedure consisting in the generalization of each new object occurring at a definite 
order in the expansion to a whole family of more complicated objects and the symbolic 
evaluation of all the members of the family insures a considerable reduction in the number 
of new objects appearing at each further step in the extension of the series. 

A final comment concerns the opportunity of applying the above strategy directly to the 
computation of connected group-theoretical factors, i.e. of the potentials. The generation 
of disconnections can be performed algorithmically; however we did not develop a specific 
computer program, resorting to geometric arguments in the cases we analyzed explicitly. All 
these cases were simple enough for us to be able to write down compact symbolic expressions 
referring directly to the potentials. Some of our results will be presented in detail in the 
following Sections. 



16 



VI. TECHNICAL REMARKS ON GROUP INTEGRATION 

In evaluating quantities like S^^y , one may take advantage of the invariance properties 
of the Haar measure for group integration 

di^{U) = di^{UA) = di^{AU) (36) 

in order to ehminate ( "gauge" ) one of the variables (defined on the nodes of the diagram) . A 
judicious use of gauging can induce notable simplifications in the actual computations, by re- 
placing "open indices" (representations) with "closed indices" (characters) in the integrands, 
and decoupling many variables from each other. 

As an illustrative example, let us consider the simplest nontrivial superskeleton. In 
principle we must evaluate 

Ck.3,r-4,.5,.6 « / X(n)(^5^)XM(5C't)x(.3)(CA^)X(M(^^^)X(r-5)(5^^)X{.6)(C^^) 

dAdBdCdD. (37) 

However, by gauging the variable D we can reduce the previous expression to the factorized 
integral 

'^n,T-2,r-3,r-4,r-B,r-6 ^ j X(r-4)(^^) 

X / Xir,){B') &t){B) ^(%(i?) dB j Xire){C^) ^^^{C) ^[l^iC) dC , (38) 

whose factors in turn will be expressible in terms of the representations of the identity via 
the relationship 

/ X(*)(^^)^J?(^)^J)(^)rf^ = E / Xit){A')^[:){A)C^rsu)Z,8dA 

= 4^s«;;,. = 4e<. (39) 

where C(rst)[^^^5 are the Clebsch-Gordan coefficients and 5^*^^^^ are the (not necessarily irre- 
ducible) representations of the identity. 

We shall call these factors "gauged vertices" , and present a few explicit examples, because 
they are essential ingredients of most of the actual computations wc have performed; it is 
immediately apparent that proper gauging can reduce the evaluation of all X (as well as Y) 
superskclctons to contractions of gauged vertices. 

The simplest nontrivial vertex involves two n = 1 representations and one n = 2 rep- 
resentations. There are three n — 2 representation, which we write down adopting the 
notation 

&'t'\A) = \[A,,Aki±AuAkj\, (40) 
^ti\A)^A,Al-^5,,5,i, (41) 
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where ^+ = ,^2;0 and = ^i,i;o- One can easily show that the (ungauged) vertices are 

{6ib6kd ± SidSkb){SajScl ± SalScj), (42) 



4d4 



J {^)^ba^cd " ^ (^ik,db"ca,jl 



The gauged vertices are trivially obtained by contraction of indices, and correspond to the 
representations of the identity matrix in the form (^OD, (|41|). Eqs. (^) and ( ^31) may also 
be used in the evaluation of a few integrals belonging to superskeletons with topology H. 

The next vertices in order of difficulty involve one each of the n = 1, n = 2, and n = 3 
representations. Adopting for n = 3 representations the notation 

(3) (3) (2,1) (2,1) f..s 

X+ = X3;0, X- = Xl,l,l;0, X+ = X2;l, X- = Xl,l;l, (44) 

we may express the corresponding vertices in the form 

/ xf{A)&t'\A^)Al^dA=^^5f{ikm,3ln), (45) 

/ X2,i;o(A) ^t'\A^) Aln dA = 6^2^'''\tkm,jln), (46) 

^(2;i)(^) ^t^^ dA = -^ 5^^'\^kn,Jlm), (47) 



J xf''\A) ^t^'-{A^) A,, dA=-^^ 5t^'\ikn,jlm), (48) 



where 



5^i\ikm,jln) = l[SijSklSmn + SuhJmj + Sinhj^ml 

± SilSkj5mn ± Sij6knSml ± ^m^fci^mj], (49) 
6±''^'°\ikm,jln) = l[26ijSklSmn - ^U^kn^mj - SinSkjSml 

± '^^ilSkjSmn T SijSknSml T ^inSklSmj], (50) 

6±'^\ikn,jlm) = I[6ij6ki6nm ± SijhmSni] 

— 2(jV ± 1) [^^^'^J'^""' ^il'^j^^nl ± ^in^jl^km + ^in^jm^kl\- (51) 

Aside from a few technicalities, the results from group integration presented in this 
section are essentially all that is needed for an evaluation of the full 15th-order strong- 
coupling contribution to the fundamental two-point Green's functions of the two-dimensional 
chiral model on the square lattice. 
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VII. COMPUTING THE POTENTIALS 



The quantities that we have denoted with the general symbol P|„ and called potentials 
are the connected parts of sums over the sets of representations consistent withy the geome- 
try of a given skeleton diagram. Needless to say, knowledge of compact analytic expressions 
for wide classes of potentials can only dramatically simplify the task of evaluating explicitly 
high orders of the character expansion. In turn, since the reduction of any diagram to its 
superskeleton can be performed algorithmically, simplifying the problem of diagram recog- 
nition, it would be obviously pleasant to possess expressions for potentials general enough 
to be referred to superskeletons instead of individual skeletons. 

We made some progress in this direction, classifying all and evaluating most of the 
skeleton diagrams whose superskeletons are drawn in Fig. ^ and obey the constraint n < 3 
for all links. In this section we shall present some general considerations and all the results 
that are needed for an explicit evaluation of all G{x) up to 12th order. We computed many 
more potentials, but often results are too cumbersome to make their presentation useful in 
any sense; they are available upon request from the authors. 

We recall that the potentials are labeled by the same symbols attributed to the corre- 
sponding skeleton diagrams. 

We already mentioned that the length of the n = 1 links does not enter the definition 
of the potentials. Moreover, the bubble content of the n = 1 links is factorized, i.e. the 
connected value of the full diagram is simply the product of the connected values of the 
diagram without bubble insertion and the diagram obtained by closing the bubble on itself 
and dividing by A^^; both these quantities are just lower-order potentials. The proof of 
factorization is very simple, and can be obtained immediately by gauging one of the vertices 
of the bubble and integrating over the second vertex variable. 

This explains why we decided not to have a notation for skeletons with bubbles along 
n = 1 paths: their name and value are expressed by the product of their factors. 

Let us now consider bubble insertions on nontrivial links n 7^ 1, in order of difficulty. 
The simplest case involves insertion of a bubble formed by two n = 1 lines between two 
vertices. Let us work out this example in detail in order to explain the general procedure. 
We take the product of representations 

((1;0)©(0;1))®((1;0)©(0;1)) 
= (2; 0) © (1, 1; 0) © (0; 2) © (0; 1, 1) © (1; 1) © (1; 1) © (0; 0) © (0; 0). (52) 

We must recognize that the existence of such a bubble imphes the possibility of two dis- 
connections of the total diagram, corresponding of the two orientations of the closed path 
running around the bubble. Therefore the connected contribution of the bubble is obtained 
by removing the two (0; 0) representations from the product, and amounts simply to replac- 
ing the bubble with a single n = 2 line (of length L = 0), with weight obtained from Eq. (p9D 
and expressible in the form 

B± = — for^i, B,,^ = -— for^i;i. (53) 

U-I- 

Given the ubiquitous presence of insertions of such bubbles along n = 2 lines, we will adopt 
the shorthand notation 
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a^2,L,b,... = 2,L, [l;l]^ 



(54) 



for the insertion of 6 [1; 1] bubbles. Such insertions imply the replacements 

d± ^ d±(S±)^ (55) 

^ rfl;l(Si;i)^ (56) 

in the expression for the value of the corresponding superskeleton, and the inclusion of a 
factor 2^ in front of all the disconnections corresponding to a splitting of the n — 2 hue. 

We may now consider the insertion of a bubble [1; 2] between two vertices. According to 
the rules for the product of representations, the corresponding contribution is obtained by 
replacing the bubble either with a n = 1 line or with a n = 3 line (with different weights 
attached to different n — 3 representations). In the first case we may apply the previously 
discussed factorization, while in the second case it is convenient to define the bubble factors 
B{a, b) as 

dfB^l\a, b) = Nd±zlBi, (57) 
^^2,i;o52,i;o(a, b) = Nd+zlBl + Nd-z''_Bt, (58) 
^(2,1)^(2,1)^^^ 6) = Nd±zlBi + Ndi.izl.^B\.^. (59) 

The insertion of the set of k bubbles [2, ai, 6i; 1]...[2, Cfe, 1] along a n = 3 line can now be 
accounted for by the following substitutions in the expression of the superskeleton: 

7(3) _ ^(3) A d(3). 



df ^dt'\lBt\a,M), (60) 

1=1 

k 

C^2,l;0 C?2,l;0 H ^2,l;o(ai, h), (61) 

i=l 

dt^'^^df'^i{B^^^'\a.M). (62) 



Moreover one must introduce factors of 2^' in the disconnections involving the splitting of 
the ith n = 2 line, and a factor 3*^ in the disconnections involving the full splitting of the 
n = 3 line into n = 1 lines. 

Next in order of difficulty are the rules concerning the insertions of [1, 3] and [2; 2] bubbles. 
In each case the allowed replacements involve either an = 2oran = 4 line. 

The bubble factors to be inserted along a n = 4 line are essentially trivial generalizations 
of our previous examples whose expressions we shall not exhibit explicitly. 

The n = 2 case is more interesting, because it is the first instance of a new phenomenon: 
the occurrence of disconnections of the skeleton diagram not corresponding to disconnec- 
tions of the superskeleton. As one may easily understand, these disconnections correspond 
to lower-order bubbles that may be removed from the skeleton turning it into another ac- 
ceptable skeleton. This possibility can be systematically taken into account by defining 
connected bubble insertions. 
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FIG. 6. Potentials. 



Let us therefore introduce the bubble factors B{p; br), corresponding to the 

insertion of [1; 3,p, [1; 2, Oi, 2, a^, br]], and C{ai, bi, a2, b2), corresponding to the inser- 

tion of [2, ai, 6i; 2, 02, 62]: 

r r 

d±B±{p- ai, bi- a,, 6,) = ^'^(^1'^ H ^± V^, + iV^^2,i;o4,i;0 H 52,i;o(ai, 

■t=i j=i 

- 2N^d±zi n (4"^± + 2^^5±) , (63) 

i=l 

di,iBi,i{p- ai, 61; a„ 6,) = Ndf'\zf'^f 

- iV^rf.^iz?^, n (^iB^ + 2''^i?i;i) , (64) 
ci±C±(ai, 61; 02, 62) = di,izl]^B\]^ Id+zl^B':^ + d^zZ'B'^ 



+ di,izllB\l 

di;iCi-i{ai,bi; a2,b2) = dl.izl^^^^"'^ B^}^^^ 



d+zl'B''^ + d^z''_'BtA - 2N^2^^^^\ (65) 



+ 



d+z^^B^^ + d-Z^'.^B^^] [d+z°^B^^ + d-Z^JBt' 
2^42bi+b2_ 



(66) 



When considering disconnections of these diagrams, one must be careful to include only 
those that have a corresponding term among the disconnections of the superskeleton. 

These rules are the essential ingredients for the construction of the connected contri- 
butions of all the skeleton diagrams entering our 15th-order calculations. In particular, all 
potentials entering 12th-order calculations can be obtained by the abovementioned insertions 
into the superskeletons drawn in Fig. M. 
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The values of these potentials are reported in Appendix ^ Here we will only report the 
results concerning W(2, ...), for reference and illustration of our formalism. 

We first recall that W(l) is completely trivial: W(l) = 1 and the associate geometrical 
factor is related to the number of self-avoiding random walks of length equal to the power 
of z. 

For the most general potential related to W(2, L) we are interested in, the main n = 2 
line splits into q [1; 1] bubbles, r [2; 2] bubbles (the bubble links themselves splitting into hi^i 
and hi 2 [1; 1] bubbles), and s [3; 1] bubbles (the n = 3 link splitting into Uj [2; 1] bubbles, 
each n = 2 link splitting into b'j,^ [1; 1] bubbles). We obtained the value of the potentials in 
the form 

A^^W(2, L, q, [2, ai,i, bix, 2, ai,2, &i,2]---[2, a^,!, Kx, 2, 0^,2, &r,2] 

X [l;3,p,,[2,a[^,XyA]---%<u,Au^-^]]-A^-,^^PsA^,<^^^ 

s r 

= z^dlBl Y[ B+{pj, a^. 1, 6^. 1, aj>., &j>^) H C+iai^i, h^i, 0^,2, &i,2) 
j=i i=i 

s r 

+ z^(f_B'L Y[ B^{pj, a^. 1, 6^. 1, a^.„,, 6^.„^) C_(ai,i, 6^,1, 0^,2, &i,2) 
j=i i=i 

s r 

+ \zi;idl,iBli n ^i;i(Pi> &j>,) n C'i;i(«i,i' «i,2, ^i,2) 

i=l i=l 
- n Y: 2^-^-^(".)''-'"W(2,p, + + Sfce^(..)(^;> - 1)) 

X n(2^-^W(2; a,,i, 6,,i) + 2''-iW(2; a,,2, ^^,2)), (67) 

i=l 

where ^{uj) are all the subsets of and Zlm052(Mj) is a shorthand for 

Sme{l,...,nj}\5^(Mj)- 



VIII. THE TWO-POINT GREEN'S FUNCTIONS AND THE INVERSE 

PROPAGATOR 

The techniques and results presented in the previous Sections set up the stage for the 
evaluation of the strong-coupling series for the two-point Green's functions G{x) of U(A^) x 
U(iV) principal chiral models on a two-dimensional square lattice, as functions of x, z{(3), 
and of the potentials. At any finite order q of the strong-coupling expansion, only a finite 
number of coordinate space Green's functions are nonzero, owing to the fact that the leading 
contribution comes from the shortest walk connecting x with the origin, which is proportional 
to 2l^il+l2;2|. therefore all the Green's functions such that |xi| -|- 1x21 > q vanish. The number 
of nontrivial Green's functions, exploiting discrete symmetries, is therefore 

\{q + 2f (geven), i(g + l)(g + 3) (g odd). (68) 

Coordinate space Green's functions are the natural output of a strong-coupling com- 
putation. Listing their individual strong-coupling series is however by no means the most 
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compact and physically most appealing way of presenting the results. It is certainly conve- 
nient to introduce the lattice momentum transform 

G(j)) = G{x) exp{ip ■ x) (69) 

X 

which, because of the lattice symmetries, turns out to be a function of the symmetric 
combinations of cosriipi and cosn2P2i with ni,?7,2 < q. 

A really dramatic simplification however occurs only when we take into consideration the 
inverse lattice propagator G'^ip). Indeed, due to the recursive nature of the path-generating 
process, any strong-coupling expansion admitting a reinterpretation as a summation over 
paths can be seen, at any definite order in the expansion, as originated by a generalized 
Gaussian model in which the appearance of new structures violating lower-order recursion 
equations can be seen as the effect of quasi-local interactions that appear in the inverse 
propagator as Fourier transforms of non-nearest neighbor couplings. A new structure capable 
of violating the recursion must correspond to a nontrivial path topology, with the property 
of multiply connecting the endpoints. Such a path must necessarily be at least three times 
as long as the minimal path. This arguments shows that, in contrast with Eq. (|69|) , in the 
inverse propagator combinations of cosnipi and cosn2P2 may appear only for ni,n2 < g/3. 
A more refined analysis shows that the highest values of rii and n2 generated in G~^{p) to 
order q in the expansion are 

ni,n2<u — 2 (g = 3M — 2), ni,n2<u — l {q = 3u — l), 

ni,n2<u (g = 3m), (u integer). (70) 

A more immediate physical interpretation of the results is obtained by introducing the 
traditional function 

p, = 2sm^, (71) 

and expressing G~^{p) as a function of = 2(1 — cos p^). One may easily get convinced that 
the number of independent symmetric combinations of powers of entering a given order in 
the expansion of G~^{p) is equal to the number of independent effective couplings one might 
define at the same order, consistently with the abovementioned considerations. This is in 
turn related to the number of lattice sites, not related by a lattice symmetry transformation, 
such that + |x2| < u. We found that a natural basis for the parametrization of these 
independent combinations is offered by 

P'-m'-py, s + 2t<u, (72) 

where 

p° = l, p'' = Y.Pl' (^>l)- (73) 

We also found that terms with t ^ appear in G~^{p) at order g = 3(s + 2t), while terms 
with t = appear only at order q = 3s + 2; the implication of this phenomenon will be 
discussed later. 
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We will therefore make use of the parametrization 

oo u 

G'\p) = Ao + A^p' + Y: E A^,sp''m"-P'f-'^^", (74) 

u=2 s=0 
u~s even 

to present our strong-coupling results for the inverse propagator in the form of expansions 
for the coefficients Au^s- As already mentioned, the expansions of Au^s will be power series 
on z, starting with z^" when s ^ u and with z^'^^'^ when s = with coefficients that are 
polynomials in the potentials. 

The large- iV limit of the A^^^ is 

= 1 - 4 2 + 4 ^2 - 4 + 12 / - 28 + 52 2^ - 132 + 324 - 908 z^ + 2020 z^^ 

- 6284 2" + 15284^^2 _ 48940^13 ^ \\^^yiz^^ - 393132 + 0{z^^), (75a) 
= ^ + ^3 _^ 7 ^5 ^ 4 + 33 + 32 + 243 + 324 + 1819 z^^ + 2520 z^'^ 



+ 14859 z^^ + 23084 z^^ + 123883 z^^ + 0{z^^), (75b) 

^2,0 = - 6 - 8 2^ - 57 z^^ - 116 - 500 z^'^ - 1152 - 5155 - 11632 z^^ 

+ 0(2^^), (75c) 
^2,2 = -2 2^ - 4 - 24 - 70 - 242 - 816 z^^ - 2824 - 8528 z^^ 

+ O(^i^), (75d) 

^3,1 = ^9 + f + 26 + 144 ^13 ^ ^14 ^ ^gQg ^15 ^ 0(^1^), (75e) 

^3,3 = 2 + 4 + 40 ^ ^40 ^14 ^ 548 ^15 ^ 0{z^^), (75f ) 

^4,0 = -| z^^ - 37^1^ - 84^15 + 0(^^'^), (75g) 

^4,2 = -2^^ - 31 - 64 z^^ + 0(^1^), (75h) 

^4,4 = -2^1^ -4^^5 + 0(^1^), (75i) 

^5,1 = 7^15 + 0(^16), (75j) 

^,3 = ^'' + 0(zi^). (75k) 



We have computed all A^^s to 0(2^^) as functions of the potentials, but the results will 
not be presented here, for reasons explained in the Introduction. We shall limit ourselves to 
the presentation in Appendix y of 15th-order expressions as explicit functions of z and N . 
These functions are obtained by substituting Eq. ([T0|) for the character expansion coefficients, 
obtaining the values of the potentials as A^-dependent coefficients, and summing up all 
homogeneous contributions. 

The limitations of such a procedure can easily be identified: gth-order expressions are 
correct for U(A^) groups with > g/2 and for SU(A^) groups with > g + 2. Even 
symbolic expressions for potentials suffer from some limitations, essentially because for small 
N not all the representations formally introduced are really nontrivial or independent. A 
manifestation of this fact is the appearance of the so-called 't Hooft-DeWit poles, which 
plague U(A^) strong-coupling expressions when < (g — 2)/4. In SU(A^) another limitation 
comes from the occurrence of self-dual representations, which spoil the applicability of U(A^) 
results already for A^ < g/2. 

In practice the results we have presented hold as they stand for all U(A^) groups with 
A^ > 7, while by using 12th-order expressions in terms of potentials one might obtain with 
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a minor effort 15th-order expressions correct for all > 3. SU(A^) groups are correctly 
reproduced for N > 16. and by use of 8th-order potentials one might obtain all N > 7. 

We must however stress that in their most abstract formulation, i.e. when expressed as 
weighted combinations of connected group-theoretical factors, our results are fully general 
and apply not only to principal chiral models but also to all nonlinear sigma models on 
group manifolds admitting a character expansion, including O(A^) and CP^~^ models. 

A number of physically interesting quantities can be extracted from G~^{p) by appropri- 
ate manipulations. In the present Section we will only present the large- limit of some of 
them. In particular we obtain the magnetic susceptibility 

^ 1 + 4:Z + 12z^ + 36z^ + 100 + 284 z^ + 796 + 2276 z'^ + 6444 / + 18572 z^ 

+ 53292 z^^ + 155500 z^^ + 451516 z^^ + 1330796 z^^ + 3904908 z^^ + 11617356 z^^ 
+ 0{z^^). (76) 

By defining the second moment of the correlation functions 

x{x')^-\i:^'Gix)^^,, (77) 
we can introduce the second-moment definition of the correlation length 

= - - 4 + 3 z + 2 - 4 + 12 - 40 / + 84 / - 296 + 550 - 1904 z^° 

z 

+ 3316 z^^ - 15248 z^'^ + 27756 z^^ + 0{z^^) (78) 
and the corresponding wavefunction renormalization 

= ^ = ^ -t- + 7/ + 4 / + 33/ + 32/ + 243/ + 324/° 

+ 1819 /^ + 2520 z^^ + 14859 z^^ + 23084 z^"^ + 123883 z^^ + 0{z^^). (79) 

The true mass gap should in principle be extracted from the long-distance behavior of 
the two-point Green's function: 

11^- lim In ^ . (80) 

|a;|^oo \x\ 

This quantity is however related by standard analyticity properties to the imaginary mo- 
mentum pole singularity of G{p), and we can therefore extract the mass gap by solving the 
equation 

G-\pi^ifXs,P2=0)^0. (81) 
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In absence of strict rotation invariance, this quantity is to be interpreted as the wall-wall 
correlation length. Let us notice that Eq. (|8l|) involves only the coefficients Au^u in the 
expansion ( |7iD of G~^{p). the power series of these coefficient in turn start with ^3ii+2 ^^^^ 



are associated with factors p^". Eq. ( pi]) is therefore, order by order in the strong-coupling 
expansion, algebraic and series-expandable in the variable 

z = 2z{cosh fis -1). (82) 

By recalling the properties of yl„ „ one can easily get convinced that knowledge of G^^{p) 
to 0(2;^""^) and 0(2;^'^) allows for the determination of Hs to 0(2^") and ©(z^""*"^) respec- 
tively. There is a deep connection between the orders of the strong-coupling expansion "lost" 
in the evaluation of fig and the above-mentioned considerations on the appearance of struc- 
tures violating the recursive relationships among paths. Indeed these structures break down 
the exponentiation of the wall- wall correlation functions at short distances [|16],|1^0, and 
one is easily convinced that loss of exponentiation at a distance ~ u implies from Eq. ( |5DD 
a residual precision ~ 2u in the determination of /i. 

The resulting value is 

/i, = -logZ-2^-|z3_2^4-f ^5_g^6_310^7_7Q^8_ 3^^9_520;2l0- 

_ 13154^12 + 0(^^3)^ (83) 

In full analogy with the discussion above, we may consider the equation for the diagonal 
mass gap (i.e. the diagonal wall-wall correlation length) 



G-\pi=tfXd/V2,p2=tfXd/V2) = 0. (84) 
Eq. ( P^ is algebraic and series-expandable in 

zM^ = Az l^cosh ^ - 1) • (85) 

Moreover one may show that zM"^ is an even function of z, and knowledge of G^^{p) to 
0(2;^") allows for the determination of to 0(z^"). The result is 

= - log 2;^ - ^2 _ 3 ^4 _ 119 ^6 _ i3g ^8 _ 4|p ^10 + Q^^12y (gg) 

Our results for the side and diagonal mass gap in terms of potentials, up to 0(2;^^) and 
0(2;^^) respectively, are available as expained in the Introduction. They are presented in 
form of explicit functions of and z in Appendix |C[ In order to compute the 12th-order 
contribution to /i^, we evaluated a few long-distance Green's functions to 0{z^^) and 0{z^'^). 

The analytic properties of the strong-coupling series (radius of convergence, zeroes of 
partition function, critical behavior) are best studied by considering such bulk quantities as 
the free energy 

F(/5) = ^In/ \{dU{x)eM-SL). (87) 

the internal energy 
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and the specific heat 



E{P) ^ 1 - ~ = 1 - G{(1,0)), (88) 



We were able to generate 20th-order series for the free energy per site F in terms of the 
potentials introduced in Fig. the results are available as explained in the Introduction. 
We shall only present here the explicit expression for large N up to 18th order, and report 
the result in terms of and (3 in Appendix p, with the usual warning that they hold only 
for V{N), N>9, and SV{N), N > 18. 

F = 2z^ + 2z^ + Az' + 19z' + 9Qz''> + 604 z'^ + 4036 z'^ + z'^ + ^ z'' 

+ 0(^2°). (90) 



For large but finite A^, we can use Eq. (^2]) to replace the derivative with respect to /3 
in Eq. (^) with a derivative with respect to z, thus obtaining a relationship between F 
and G{{1,0)), which we verified explicitly. It should be noticed that no simple relationship 
between F and G{{1, 0)) can be obtained in terms of potentials. 



IX. STRONG COUPLING EXPANSION ON THE HONEYCOMB LATTICE. 

On the honeycomb lattice we consider the action with nearest-neighbor interaction, which 
can be written as a sum over all links of the lattice 

5h = -2N(3 Tr [Ui 0] , (91) 

links 

where /, r indicate the sites at the ends of each link. As on the square lattice, the link length 
a is choosen as lattice spacing, i.e. as length unit. The continuum action of chiral models is 
obtained from the a —>■ limit of by identifying 

T = ^ . (92) 

For what concerns the strong-coupling expansion, we only mention that the determina- 
tion of the geometrical factor is straightforward in the light of our general discussion. The 
only difference with the square lattice concerns the generation of non-backtracking random 
walks. It is easy to see that the honeycomb lattice can be mapped in a subset of the square 
lattice, having the same sites, the same links in y direction, and only the links in x direc- 
tions starting from even sites: therefore the walks on the honeycomb lattice are a subset 
of the walks on the square lattice. The value of the potentials is obviously unchanged, but 
since the computation can be pushed some orders further on the honeycomb lattice a few 
new calculations are needed. The only subtle point is that, since only half of the sites of 
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the lattice are related by translation invariance, the free energy per site is one half of the 



quantity computed according to Sect. |IV. 



We generated strong coupling series of the free energy up to 0{z'^^), and of the fundamen- 
tal Green's function up to 0(2^°). Our results as functions of the potentials are available 
as explained in the Introduction. We present here only the large-iV results; we refer to 
Appendix ^ for results for large but finite A^. 

In analogy wityh the square lattice, we evaluated the strong coupling series of the fun- 
damental correlation function G{x) = (Tt [U{x)U^ {0)]^ as function of x, z{f3), and of the 
potentials. The number of nontrivial components of G{x) which must evaluated at a given 
order q is 

if g = 4A;,4A; + l, 

^ ^ ^ if g = 4A; + 2,4A; + 3, (93) 

with k integer. 

The magnetic susceptibility x ^-nd second-moment correlation length E,g are defined on 
the honeycomb lattice in perfect analogy with square lattice definitions: x = J2xG{x), 
xeG = jE.x'G{x). 

The analysis of models on honeycomb lattices presents some complications, which will be 
illustrated in some detail in Appendix |^ by considering a simple Gaussian model of random 
walk. The point is that, unlike square and triangular lattices, not all sites are related by 
a translation; this fact does not allow a straightforward definition of a Fourier transform. 
Only sites at an even distance (in the number of links) are related by a translation. We 
therefore define even and odd fields Ue, Uo] Ue{x) = U{x) for even x, zero for odd x, and 
Uo{x) = U{x) for odd x, zero for even x (the parity is defined with respect to an arbitrarily 
chosen origin). We then define and even and odd correlation functions 

Ge{x -y) = l^ (Tr [U,ix)Ul{y)]) = {Tt [f/o(x)f/t(y)]) , 

Go{x -y) = j^ (Tr mx)UKy)]) = 1 (Tr [f/o(x)f/t(y)]) • (94) 

Since even and odd sites lie on two distinct triangular sublattices, it is possible to define 
consistent Fourier transforms on each sublattice. 

Guided by the analysis of the Gaussian model, we considered two orthogonal wall-wall 
correlation functions: 

G^;^\x) = Y.G,{x,y), (95) 

y 

G^'\x) = E [G,{x,y) + Go{x,y)] . (96) 
y 

In the strong-coupling domain both G^™^ (x) and G'^'^ (x) enjoy exponentiation for sufficiently 
large lattice distance, allowing the definition of two corresponding masses /ii and fi2- 
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G'i^^(x) (xexp(-|//ix), (97) 
G^ ix) oc exp(— \V?>H2x). (98) 

In the continuum limit jii = ^2 and they both should reproduce the physical mass M propa- 
gating in the fundamental channel. The ratio Hi/ H2 allows a test of rotational invariance, in 
analogy with the side/diagonal mass ratio of the square lattice. It is also possible to define 
the quantities 

zMl = |^(cosh - 1), (99) 
zM^ = |^(cosh i\/3/X2 - 2), (100) 

which play the role of zM^ and zM^ in the determination of imaginary momentum pole of 
the inverse Fourier-transformed Green's function. 

In full analogy with the square lattice, we define the magnetic susceptibility 



X = $^G(x) = 1 + 3z + 6^^ + 12^^ + 24^^ + 48^^ + 90^^ + 174^^ + 348 + 702 



z^ 



+ 1392 + 2814 + 5658^^2 + 11532^^^ + 23706 + 49368 + 101436^^^ 

+ 21 1290 + 440598 + 928614 + 1950390 ;22° + 0(;2^^), (101) 



the second moment of the correlation functions 

) =i 

'g 4 



X = 7 E{(9^i + ^4)Ge{x) + [(3X1 - 1)' + 3xl]Go{x)}, (102) 



X 



and the second-moment definition of the correlation length 

M^ = -l—^lz-^-A + lz-Sz'' + fz'' + 8z^-16z^- 88z^° + 96z^^ 

- 144 z^^ + ^z^^ - 40 z^^ - 200 z^^ - 5520 z^^ + 5848 z^^ - 4208 z^^ + 0{z^^). (103) 

We were able to obtain 26th-order results in terms of potentials for the free energy per 
site: 



F^lz' + z' + Zz'^' + Qz'^ + 12 z^^ + 114 + ^ z^' + 1080 + 5754 z^^ 

The internal energy (per link) and the specific heat can be obtained by 

^(/5) = 1-^^ = 1-G((1,0)), (105) 

The same caveats of the square lattice case apply to the relationship between F and G((l, 0)). 
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APPENDIX A: VALUES OF SELECTED POTENTIALS 



In Sect. |Vll| we reported a rather general form for potentials related to W(2, L) by bubble 
insertions. Here we present all potentials needed for a 12th-order computation, and some 
generalizations. We also list all other potentials we computed; their expressions are too long 
and cumbersome to be reported here, and they are available upon request from the authors. 

Let us introduce a shorthand notation for bubble insertion. A sequence of r [1; 2] bubbles 
will be denoted by 

yriai,bi,...,ar,br) = [l;2,ai,bi]...[l;2,ar,br.]; (Al) 

the arguments ( br) will often be left understood. A sequence of r [2; 2] bubbles 

and s [1; 3] bubbles (with [1; 2] splittings along the n = 3 line) will be denoted by 

<^r,s 0,2,1, b2,i, Oi^r, bi,r, 02,r, b2,r]Pi, o\,b\, , ; . . . ; p<j, , 6^ , . . . , , 6^ J = 

[2, ai_i, 2, 02,1, &2,l]---[2, ai,r, &l,r; 2, 02,r, &2,r] 

X [l;3,pi,A.(ai,6i,...,ai^,6;,J]...[l;3,p„A.K,6^...,<^,6:j]; (A2) 

the arguments (ai^i, ■■■,b^J will also be left understood. 
The shorthand notation 

=4)^(0- (A3) 

will also be used. 

Potentials related to W(3, L) are easily obtained from the identity 

W(3, L, J^r) = W(2, Or, br, [1; 3, L, J^r-i]). (A4) 

One must however compute explicitly the r = case 

W(3, L) = ^ \^i'\L, 2) + 2) + ^2,i;o{L, 2) + 2) + 2] 

-2N^W{2,L)--N^. (A5) 
3 

We have also computed W(4, L, J^r,s), whose structure is similar to W(2, L, b, J'r,s) pre- 
sented in Sect. |V11| . More generally, we must mention that a generating functional for 
W(n, V) can easily be constructed by exploiting the general strong-coupling solution of the 
chiral chain problem [|T6],|T^0. It is easy to get convinced that 

^ In E E 2)^"^ = E W(n, V)z-\ (A6) 

where the sum on the l.h.s. is extended to all representations with the same value of n. 
In principle, generating functionals for potentials with bubble insertions can also be con- 
structed. 

Let us now consider the first nontrivial superskeleton Y: 
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Y(2, ai, 61; 1; 1; 1; 1; 2, as, h) = z^}+^^ B'^^+'' + ^N^ z'^^.B'l^^^ [zI^B'^+^ + z^JB'_^+\ 

+ iN^z'i^^Bl^^^ (415^1+1 + 2^15^1+^) - 2N^2^'+^''- 



a\ ,a2 • 



(A7) 
= 0, 

We have also computed more general objects of the form 



this quantity was first introduced in Ref. for the special choice hi 
and it was termed W, 



Y(2, ai, 61, J^r,,s^] 1; 1; 1; 1; 2, 03, &2, 



Y(l; 1; 1; 2, ai, bi, 2, 02, 62; 2, 03, 63) 
= J^i;i(ai + a2 + as, - 

+ {^^^\5,^^i[(^+(a2,l)i?t^ + Sl(a2,l)S^^)(^+(a3,l)5t' + ^„(a3,l)S^^) 
- (ir+(a2 + a3, 1)5^^+'^ + ir_(a2 + 03, 1)5^^+'^)" 
- 2N^2''^+^m{2, ai, 61)} + permutations of indices - 4 A^^ 2 61+62+63. 

was termed Wai,a2,a3 in Ref. 0. 



(A8) 



the case bi = h 



Y(3,p;2,ai,6i; 1; 2, 02, 62; 1; 1) 



N 



2fl^\p, l)zl'+''^Bl'+^^ + 2fl^\p, 
+ i^2,i;ob, l){zl'^^'B^_^+^^ + 3zl'B^^z^JB''_^ + Szl^B'l'zl'Bt' + 5^1+^2) 
+ ^l^''^\p, + z'^^B'l^izl'Bl' + 

+ ^l^''^\p, l){z^liB^i^iz'^^Bt' + z'^iB\lz''_} B^2 + 

- 2A^2W(2, p + ai, 61)2^2 - 2A^2w(2, p + 02, 62)2^^ 

- 2Ar2W(2, p, 0)2''i+^2 - 4A^426i+62 _ (A9) 

We also computed a few generalizations: Y(3,p, J^,,; 2, ai, 61; 1; 2, 02, 62; 1; 1), 
Y(3,p, J^r;2,ai,6i;2,a2,62; 1; 1; 2, 03, 63). 

Finally, we computed, for a few special values of the indices, Y(3,pi; 2, ai; 1; 2, 02; 1; 3,p2)- 

The second nontrivial superskeleton is X. We computed 

X(2,ai,6i;l;2,a2,&2;i;i;i;i;i) 

X (j^+(a2, 1)5^^ + J^_(a2, 1)5!„^ + ^i;i(«2, 

/Y^ 4 
+ — rfi;i(j^,(ai, -1)5^^ + J^_(ai, -l)B'^ + — J^i;i(ai, 



X (^+(a2, -1)5^^ + ir_(a2, -l)B'_^ + ^ ^i;i(a2, -1)5^^1 
+ N^f^i,i{ai + 02, + 4) 



- 2A^2[2^^W(2,a2,&2) + 2'2W(2,ai,6i)] -2''^+''2(4Ar* + 2Ar'). 
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61+62 / 



(AlO) 
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APPENDIX B: LIST OF POTENTIALS 



We list here all the potentials appearing for the first time at each order of the strong- 
coupling expansion; we identify potentials differing only for the values of L and for the 
number b of [1; 1] bubbles provided that 6 > (i.e. we don't identify 6 = with 6 7^ 0). In 
the following formulae, L will indicate a generic value, and multiple occurrences of L in the 
same expressions indicate any combinations of Ls (i.e. they can take different values). 

5th order: 

W(2,L,6). (Bl) 

9th order: 

W(3, L, [1; 2, L]), W(3, L, [1; 2, L, b]), Y(2, L; 1; 1; 1; 1; 2, L). (B2) 
10th order: 

Y(2,L;2,L;1;1;2,L;1). (B3) 

11th order: 

W(2, L, [2, L- 2, L,b]), Y(3, L; 2, L; 1; 2, L; 1; 1), Y(2, L; 1; 1; 1; 1; 2, L, b). (B4) 
12th order: 

W(2, L, h, [2, L- 2, L]), W(3, L, [1; 2, L]'), Y(2, L; 2, L; 1; 1; 2, L, 6; 1), 

Y(2, L, b- 1; 1; 1; 1; 2, L), X(2, L; 1; 2, L; 1; 1; 1; 1; 1). (B5) 

13th order: 

W(3, L, [1; 2, L, [1; 3, L]]), W(4, L, [1; 3, L]), W(4, L, [2, L; 2, L, 6]), W(2, L, 6, [2, L; 2, L, &]), 
Y(3, L; 2, L; 1; 2, L, 6; 1; 1), Y(3, L; 2, L; 2, L; 1; 1; 2, L), Y(2, L; 1; 1; 1; 1; [1; 3, L]), 
Y(2, L, 6; 1; 1; 1; 1; 2, L, 6), Y(2, L, 6; 2, L; 1; 1; 2, L; 1), X(2, L; 1; 1; 1; 1; 1; 2, L; 2, L), 
H(2,L;1;2,L;1;1;1;1;1;2,L). (B6) 

14th order: 

W(3, L, [1; 2, L, 6, [1; 3, L]]), W(4, L, [1; 3, L, [1; 2, L]]), W(2, L, [1; [1; 2, L]']), 
W(2, L, [2, L, &; 2, L, &]), W(3, L, [1; 2, L] [1; 2, L, 6]), Y(3, L; 3, L; 2, L; 1; 2, L; 1), 
Y(3, L; 2, L, 6; 1; 2, L; 1; 1), Y(2, L; 1; 1; 1; [1; 2, L]; 2, L), Y(2, L; 2, L; 1; 1; [1; 3, L]; 1), 
Y(2, L; 2, L, fe; 1; 1; 2, L; 1), Y(2, L, [1; 3, L]; 1; 1; 1; 1; 2, L), X(l; 1; 1; 1; 2, L; 2, L; 2, L; 2, L), 
X(2, L; 1; 2, L; 1; 3, L; 1; 1; 1), X(2, L; 1; 2, L, 6; 1; 1; 1; 1; 1), X(2, L, &; 1; 1; 1; 1; 1; 2, L; 2, L), 
H(l; 2, L; 1; 2, L; 1; 1; 1; 1; 2, L), H(2, L; 2, L; 1; 1; 1; 2, L; 1; 2, L; 1), 

R(1;1;2,L;1;1;1;2,L;1;1). (B7) 
15th order: 



32 



W(3, L, [1; 2, L, [2, L; 2, L]]), W(3, L, [1; 2, L, [2, L; 2, L, 6]]), W(3, L, [3, L; 2, L, 6]), 

W(3, L, [[1; 2, L]; 2, L]), W(4, L, [2, L, 6; 2, L, 6]), W(3, L, [1; 2, L]^), 

Y(3, L; 2, L; 1; [1; 3, L]; 1; 1), Y(3, L; 2, L; 1; 2, L; 1; 3, L), Y(3, L; 2, L; 2, L; 1; 1; 2, L, 6), 

Y(3, L; 2, L, 6; 1; 2, L, 6; 1; 1), Y(3, L; 2, L, 6; 2, L; 1; 1; 2, L), Y(4, L; 3, L; 2, L; 2, L; 1; 1), 

Y(2, L; 1; 1; 1; 1; [2, L; 2, L]), Y(2, L; 1; 1; 1; 1; 2, L, [1; 3, L]), Y(2, L; 2, L; 1; [1; 2, L]; 2, L; 1), 

Y(2, L, b- 1; 1; 1; 1; [1; 3, L]), Y(2, L, b; 2, L; 1; 1; 2, L, 6; 1), Y(3, L, [1; 2, L]; 2, L; 1; 2, L; 1; 1), 

Y(2, L, [1; 3, L]; 2, L; 1; 1; 2, L; 1), X(l; 1; 1; 1; 2, L; 2, L; 2, L; 2, L, b), 

X(3, L; 2, L; 1; 2, L; 1; 1; 1; 1), X(2, L; 1; 1; 1; 1; 1; 2, L; 2, L, 6), 

X(2, L; 1; 1; 1; 3, L; 1; 2, L; 2, L), X(2, L; 2, L; 1; 1; 1; 2, L; 1; 2, L), 

H(3, L; 2, L; 1; 1; 2, L; 1; 1; 2, L; 1), H(2, L; 1; 2, L; 1; 1; 1; 1; 1; 2, L, 6), 

H(2, L; 1; 2, L; 3, L; 1; 1; 1; 1; 2, L), H(2, L; 1; 2, L, 6; 1; 1; 1; 1; 1; 2, L), 

H(2, L; 3, L; 2, L; 1; 1; 1; 1; 1; 2, L), H(2, L; 2, L; 2, L; 1; 1; 2, L; 2, L; 1; 1), 

L(l; 2, L; 1; 2, L; 1; 2, L; 1; 1; 1; 1), L(2, L; 1; 1; 1; 2, L; 1; 1; 1; 1; 1), 

L(2,L;1;2,L;1;1;2,L;1;1;1;1). (B8) 



APPENDIX C: SQUARE LATTICE RESULTS FOR FINITE N 

We list in the present Appendix the values of the quantities defined in Sect. [VIII] without 
further comments. The definition of these quantities and the range of validity of the results 
presented here is discussed in Sect. ]VIII. 



7-13 + 22iV2-llA^4 ^ 101 - 158A^2 ^81A^^ 

+ 12 z^ 5 +4z^ 5 

(Ar2 - 1)2 (Ar2 _ i)^ 

q -243 + 358A^2_ 227A^^ 
+ 4z 



+ 4z 



10 



{N^ - if 

2444 - 6391 + 7193 A^^ - 3577 A^^ + 505 A^^ 

\3 



{m - 4) {N^ - ly 

+ 4 z^^ (-23568 + 89928 A^^ - 154753 A^^ + 145700 A^^ 
- 74058 A^^ + 17674 N^^ - 1571 A^^^) 
1 



{m - 4Y {N^ - ly 

+ 4 z^^ (59312 - 227352 A^^ + 401499 A^^ - 379396 A^^ 
+ 189036 A^^ - 43842 N^^ + 3821 A^^^) 
1 

^ (Ar2 _ 4)2 _ 1)4 

+ 4 z^^ (-145008 + 560056 A^^ - 1018759 A^^ + 986360 A^^ 
- 519366 A^® + 128954 A^^° - 12235 A^^^) 
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1 

+ 4 z^^ (3296016 - 16590056 + 38621009 A^^ 

- 50952462 A^^ + 40255375 - 18791194 7V^° + 4822317 A^^^ 

- 617438 7V^^ + 29153 7V^^) 

1 



X 



(Ar2 - 9) (iV2 - AY {m - ly 

+ 4 z^^^ (-73339344 + 453776040 A^^ - 1300249437 A^^ 

+ 2195554892 A^*^ - 2367482622 A^^ + 1669383326 A^^° - 759240644 A^^^ 
+ 212701682 A^i^ - 34650066 A^^^ + 2929064 A^^^ - 98283 A^^") 

X 9 9 f,+0(z^^), (CI) 



3 ^5 .6^' + ! o 7 13-22Ar2 + llA^4 

o -5 - 2A'- + lA^-^ „ 251 - 291A^2 ^ 243A^^ 
+ 8z^ 5 h z^ 5 

(Ar2 _ 1)2 (7V2 _ 1)2 

in -284 - 13 Ar2 + 441 A^^ - 369 A^6 + 81 A^s 

+ 4:Z 

{N^ - 4) {N^ - if 
+ z^^ (24208 - 76936 + 129833 A^^ - 134180 A^^ + 76114 A^^ 
- 19562 AT^o + 1819 Ari2) 
1 



(Ar2 - Af (A^2 _ 1)4 

+ Sz' 



12 -3632 + 2616 A^2 + 3555 A^^ _ 3734 jye ^ 7953 _ 2821 iV^o + 315 N^^ 



(Ar2 - Ay {m - ly 

+ (148848 - 456760 A^^ + 879351 A^^ - 922184 A^^ 

+ 540230 A^^ - 146906 + 14859 A^^^) 
1 



X 



(A^2 _ 4)2 (^2 _ 1)4 



+ 4 z^"^ (-422352 + 796520 Ar2 - 475949 A^^ - 864192 A^^ 

+ 2055327 A^^ - 1703150 N^^ + 634345 A^^^ - 102496 A^^^ + 5771 A^^^) 
1 



X 



(A^2 _ 9) (iV2 _ Af - if 
+ z^^ (75941712 - 403857000 A^^ + 1114656565 A^^ 

- 1899839324 A^^ + 2117909054 A^® - 1584540590 A^^° + 779561548 A^^^ 

- 235498610 A^^^ + 40734002 A^^^ - 3591112 A^^^ + 123883 N^°) 

^ +0(z'% (C2) 



(A^2 _ 9)^ (7V2 - Ay (7V2 - 1) 
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CI — ~z — 1-2 2; T 4:Z T 

10 212 - 25 iV2 - 355 iV^ + 255 - 57 A^^ 
^ ^ {N^ - 4) (iV2 - 1)3 

, 11 -5-100iV2 + 137iV^-12 7V6-29 7V8 

12 2608 - 1944 + 195 + 1992 - 4278 iV^ + 2001 - 250 A^^^ 

- Af - 1)^ 
60 + 1485 - 1735 A^^ + 520 + 396 - 144 N^^ 



+ 2z 



+ — —^^ -4 



- 4) {N^ - ly 

+ z^^ (303120 - 666632 A^^ + 1554257 A^^ - 1748152 A^^ 

+ 205601 A^^ + 831342 A^^^ - 468499 A^^^ + 86566 A^^^ - 5155 A^^^) 
1 



^ (A^2 _ 9) (7V2 _ (Ar2 - \f 
+ 4 z^^ (-3984 - 68344 A^^ + 247735 A^^ - 387636 A^^ 

+ 321005 A^^ - 97944 N^^ - 30420 A^^^ + 21362 A^^^ - 2908 A^^^) 



{N^ - Ay {N^ - 1) 
_ .8^^2 + 1 1 + 87V2 + 27V^ 

-9-7Ar2 + 9A^^- 12A^6 

+ 22;^° 



+ 2z 



11 



(Ar2 - if 
-5 - 112 A^2 + 155 _ 12 A^e _ 35 



{N^ - if 
12 64 + 54 Ar2 + 100 N'^ - 52 - 121 
+ (Ar2 - 1)^ 

^ 13 60 + 1917 7V2 - 1843 A^^ + 640 A^^ + 606 A^^ - 204 A^^o 
^ ^ (Ar2 - 4) (Ar2 - 1)^ 

+ 2 (-7392 + 8144 A^^ _ 55774 + 92342 - 35882 A^^ 

- 16810 iV^° + 10745 N^^ - 1412 AT^^) 
1 

X 



(Ar2 - Af {m - If 
+ 2 z^^ (-4336 - 96968 A^^ + 298881 A^^ - 461772 A^^ 

+ 413767 A^^ - 131684 A^^° - 48836 A^^^ + 32296 A^^^ - 4264 A^^^) 

X T fi + 0(z^^), 

(A^2 _ 4)2 (^2 _ 1)6 V h 



9l + 8Ar2 + 2A^^ ii5 + 100Ar2_i37^4^12 7v6 + 29A^^ 
' 2(Ar2_i)2 2 - if 
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^ IP UN^ + 93N^ + 13N^ 
+ 2z^^ — 



(iV2 - 1)3 

-60 - 1485 + 1735 - 520 - 396 + 144 N^^ 

14 11 - 25 + 1585 - 2366 iV^ + 714 + 241 TV^o 

(Ar2 - if 

+ (4368 + 68152 A^^ - 230047 iV^ + 463796 - 576029 A^^ 
+ 309400 - 21332 N^^ - 20786 N^^ + 3612 A^iS) 
1 

^ 2 (Ar2 - 4)^ {N^ - if 



I3 3 = 2 z^^ -IT + 4 z^^ 

(A^2 _ 1)2 _ 1)3 

1, 18 ^ 5 A^^ + 10 A^^ 
+ 4^13 _ 



- 1) 

. . 1 + 5 

+ 4z 



14 1 + 5 Ar2 + 223 A^4 - 334 A^^ + 106 A^^ + 35 A^i° 



15 5 + 300 - 255 A^^ + 866 A^^ - 1612 + 455 A^^° + 274 A^^^ 



(Ar2 - 1) 

5 + 300 A^2_255/ 
+ 0{z''), 



^^4:N^ + 33 + 5 
^4,0 - -z — — —3 



+ z 



2 (A^2 _ 1) 

_0 _L 1 

14 



-2 + 10 Ar2 - 241 N'^ + 362 A^^ - 105 A^^ - 37 A^^o 



(Ar2 - 1) 



5 



1, l + 2A^2^39Ar4 + 296A^6 + 28A^s 
3z^^ 1 ^0(2;^^), 



M,2 ^ -Z 



12 2 Ar2 + 9 Ar4 + A^6 



(Ar2 - 1)^ 



14 -l-5Ar2- 207Ar4 + 306Ar6-98iV8-31iV^° 

+ Z E 

(Ar2 - if 
7V2 _ if ^ 



(A^2_l)3 (A^2_l)4 ^ ^' 
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1 _ 91 /V^ 

-log(.)-2.-i.3_2.4 + 2.^^^^^ 
6-5 + 6A^2_4^4 ^6 + 114Ar2_i55Ar4 



+ 2^° h2^ , „ 

„ 22 - 64 - 

+ 2^' 



22 - 64A^2^58 A^^ - 35 



(N^ - 1) 



3 



9 -101 - 1286 ^ 97i8A^i - 15331 A^*' + 8810 A^^ - 1591 A^^^ 
^ ^ 9 (Ar2 _ 4) (Ar2 _ 1)4 

10 -1600 + 8768 - 14596 A^^ + 13432 - 8645 A^^ + 2500 - 260 A^^^ 

+ (Ar2 _ 4)2 (^2 _ 1)4 

+ 2 (-864 - 20592 A^^ + 266170 A^^ - 740050 A^^ 

+ 869302 A^^ - 508486 A^^° + 141781 A^^^ _ 14389 A^^^) 
1 

^ 11 (Ar2 - 4)' (A^2 _ 1)5 
+ 2 z^^ (-23392 + 184400 A^^ - 571030 A^^ + 920690 A^^ 

- 927373 + 643901 - 289681 iV^^ + 68171 A^^^ _ 5577 jyie) 

X 1 f. + 0(z^^), (C12) 

3 (Ar2 - 4)' (Ar2 - 1)6 ^ ^' ^ ^ 

, , 2 42-3A^2 . -79 + 248A^2_ 238A^4 
-log (2-)--^ + -^!^ + -= ^^f^. 

8 908 - 5299 AT^ + 12080 A^^ - 13391 A^^ + 6844 - 1088 A^^o 
^ ^ 8 (Ar2 - 4) (Ar2 - 1)4 

+ z^° (172816 - 1265528 A^^ + 3877481 A^^ _ 6436255 A^^ 
+ 6193060 A^^ - 3331306 N^^ + 864838 A^^^ _ §3925 ^v^^) 

X \ . + 0(z^^), (C13) 

80 (Ar2 - 4)' (Ar2 - 1)^ ^ ^' ^ ^ 

X= 1 + 4 ^ + 12 ^2 + 36 + 100 / + 284 / 

6 -195 + 199A^2 7 543- 1106iV2 + 569A^4 



- 1 (A^2 _ 1)2 

. 49 

+ 12 



g 493- 1022 AT^^ 537 Ar4 



(A^2 _ 1)^ 



37 



, q 4067- 8550 + 4643 A^^ 



{m - 1 



i2 



in 44100 - 149869 A^^ + 182355 - 90083 A^^ + 13323 A^^ 
+ 4 2;^° 



\3 



- 4) {m - 1)' 

+ 4 (481168 - 2256648 A^^ + 4239673 A^^ - 4010564 A^^ 

+ 1961034 - 452890 N^^ + 38875 A^^^) 
1 



X 



(Ar2 - 4)' (Ar2 - 1)^ 



+ 4 z^"^ (1299728 - 6161736 A^^ + 11712625 A^^ 

- 11228076 A^'' + 5568300 A^^ - 1301614 iV^° + 112879 A^^^) 



X 



(A^2 _ 4)2 (^2 _ 1)4 



+ 4 z^^ (3526000 - 16894776 A^^ + 32535287 A^^ 

- 31650616 + 15957638 A^^ - 3784938 A^^° + 332699 A^^^) 
1 



X 



(Ar2 - {N^ - 1)^ 



+ 4 z^^ (85479984 - 509240632 A^^ + 1277936387 A^^ 

- 1742114314 A^*^ + 1385750301 A^^ - 641889782 N^^ + 163736431 A^^^ 

- 20707322 N^^ + 976227 A^^^) 

1 

X 



(Ar2 - 9) (Ar2 - Ay - ly 
+ A z^^ (2078977104 - 14813458920 A^^ + 45842748421 A^^ 

- 80400814700 A^^ + 87447748126 A^^ - 60663055822 N^^ + 26640570340 A^^^ 

- 7157578626 A^^^ + 1112718466 A^^^ - 90043336 A^^*^ + 2904339 A^^") 

X 9 9 K + 0(z^^), (C14) 

(Ar2 - 9)' (A^2 _ 4)2 (^2 _ 1)6 ^ ^ ^' y ^ 

..2 1 . o ^ s . 4^^ + l o . -5-4Ar2 + 6Ar4 

M^= z~^ -A + 3z + 2z^ -Az^ — + 2 z^ 5 

G Ar2-1 (A^2_i)2 

„ .6 + 2N^-5N^ , 7 -32 - 42 Ar2 + 21 AT^ 

+ 8^^ 5 VAz^ 5 

(iV2 - 1)2 (Ar2 _ 1)2 

8 184 + 98Ar2 - 376A^^ + 203Ar6-37Ar8 

^ ^ (Ar2 - 4) (Ar2 - if 



+ 2z^ 



-6576 - 3656 N"^ + 30629 A^^ - 32810 A^^ + 14735 - 3245 N^^ + 275 N^"^ 

(Ar2 - Af (iV2 - 1)' 



10 4016 + 4136 N"^ - 16789 A^^ + 15608 A^^ - 7937 A^^ + 2257 N'^^ - 238 A^^^ 

(Ar2 - Af {m - if 
+ 2 z^^ (-34560 - 63904 A^^ + 148416 A^^ - 142210 A^^ 
+ 74870 - 19631 Ar^° + 1658 A^^^) 
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1 

^ {N^ - Af - 1)' 

+ 8^^^ (186048 + 87904 iV^ - 727204 iV^ + 1149812 A^^ 

- 1175339 A^^ + 716157 Ar^° - 232150 TV^^ + 34766 A^^^ - 1906 7V^^) 

1 

^ {N^ - 9) (iV2 - Af - if 
+ 2 z^^ (-14774400 + 8055072 + 58712088 iV^ 

- 170052698 + 261060513 - 232073848 N^^ + 117200494 N^^ 

- 33867953 N^^ + 5360045 A^^*^ - 433255 A^^^ + 13878 A^^°) 

X 5 ^ ^ + 0(z^^), (C15) 



2 4 . „ 14-30 A^2^ 19 A^^ 



+ 8^^° 



(A^2 _ 1)^ 

7- 16A- + 12A^4 
- if 



+ 2 z^^ (5952 - 30144 A^^ + 53364 A^^ - 67746 A^^ + 37500 A^^ 

- 9589 A^^° + 906 A^^^) 
1 



X 



3 {m - Af (A^2 _ 1)4 
+ 4 z^^ (4704 - 25200 A^^ + 55366 A^^ - 62068 A^^ + 36768 A^^ 
- 10012 A^^° + 1009 7V^2) 
1 



X 



- Af {N^ - if 



+ z^^ (15230592 - 119651328 A^^ + 409792072 A^'^ 

- 799817292 A^^ + 974422200 - 760569676 A^^° + 375693595 A^^^ 

- 112460534 A^^^ + 19258826 A^^*^ - 1690814 A^^^ + 58471 N^°) 



X 



2 (N^ - 9f (N^ - Af (N^ - If 

+ 8 z^^ (237447936 - 2045848320 A^^ + 7796836128 A^^ 

- 17299894704 AT^ + 24693730379 A^^ - 23648019056 A^^° + 15403609647 A^^^ 

- 6787736700 A^^^ + 1995934103 A^^^ - 381749639 A^^^ + 45159907 A^^° 

- 2974083 N^^ + 82898 A^^^) 

X 2 4 6 + Oiz^''). (C16) 

3 (Ar2 - 9)' (Ar2 - 4)^ (Ar2 - 1)^ ^ ^ ^ ^ 
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APPENDIX D: THE GAUSSIAN MODEL ON THE HONEYCOMB LATTICE 



There are a few subtleties in the analysis of models on the honeycomb lattice, that are 
best illustrated by considering a simple Gaussian model of random walks. The essential 
point is related to the fact that lattice sites are not all related by a translation group: only 
points at an even distance (in the number of lattice links) arc related by such a symmetry. 
As a consequence, it is convenient to define even and odd fields 0e) 0oj according to the 
parity of the corresponding lattice sites with respect to an arbitrarily chosen origin, and 
even and odd correlation functions Ge = {4>e4>e) = {4>o4>o), Go = {<Pe<Po) = {4>o4>e)- 

Let us represent the Cartesian coordinates of the (finite, periodic) lattice sites by 



m, T^VSn I a, 



1, 

2 



{x, y)o = (|m + 1, \V^'n)j a. 



(Dla) 
(Dlb) 



where m and n are integer numbers satisfying the conditions 0<m<Li,0<n< 2L2, 
and m + n is even. The total number of lattice points is 2L1XL2, the number of links is 
3L1XL2, and the number of plaquettes is LixL-2. 

The finite-lattice Fourier transform is consistently defined by 



Ge{p) = E e^^"Ge(x), 



(D2) 



and similarly for G'o(p); the set of momenta is 

27r /2 m 1 n 



(D3) 



with fh and h integers, and 0<m<Li, Q < h < L2. 

In a random walk model where walks of length v are weighted by a factor , it is easy 
to establish from the recursion relations the following relationships: 



Ge{p) 

Goip) 



1 



2 cos ^VSp; 



■2 6 



3ipi/2 



l + 2cosiV3p2e-^^^i/^ 



Goip) - 

Geip), 



(D4a) 
(D4b) 



where Ge{p) and Goip) are the Fourier transform of Ge{x) and Go{x) respectively. As a 
consequence, we obtain the even momentum-space Green's functions in the form 



Geip) 



1 - /?2 1 + 4 C0S2 ^V3p2 + 4 COS IV^P2 COS I 



Pi 



(D5) 



The critical value of P is easily found to be /3c = |, and as a consequence we find the 
massless lattice (even) propagator 



Kip) 



8-4 cos^ lVSp2 + 4 cos lV3p2 cos fpi p-*o 



(D6) 



The odd propagator is simply 



40 



Ao(p) 



3*^ 



1 + 2 COS hV3p2 e~ 



■3ipi/2 



(D7) 



The structure of the propagator in the Gaussian model offers an important indication 
about the possibihty of exponentiation in wall-wall correlations. Let us indeed recall that 
exponentiation corresponds to a simple structure 

(D8) 



A — B cos p 

in the corresponding propagator. Let us now observe that Eq. ( P5|) implies 



1 - /32 5 + 4 COS f j9i 



1-/3 l + 2cosiy3p2 



(D9) 
(DIO) 



We have therefore the possibility of defining two different exponentiated "wall-wall" corre- 
lation functions, i.e. 



y 



(Dll) 
(D12) 



Even in more general models, in the strong-coupling domain, for sufficiently large lattice 

If we 



distance exponentiation will hold for the correlation functions G^^\x) and G2 
take into account the discrete rotational symmetry of the honeycomb lattice we may easily 
recognize that the above correlations can be referred to directions differing by a 7r/6 angle, 
and this is the maximal violation of the full rotational symmetry one can find on this 
kind of lattice. Therefore the ratio of the two different correlation lengths one may define 
is an optimal measurement of the violation of rotational invariance in the model under 
examination, in analogy with the side/diagonal mass ratio of the square lattice. 



y(w) 
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APPENDIX E: HONEYCOMB LATTICE RESULTS FOR FINITE N 



We list in the present Appendix the values of the quantities defined in Sect, ^without 
;nts. 1 

log 2z^ — z 



further comments. The definition of these quantities presented here is discussed in Sect. fX 



/il 



6 (A^2 _ 



9-27N'^ + 35N'^-37N^ 
4 - if 



+ z 



10 



-194 + 1136 - 2539 A^^ + 2546 A^*^ - 1174 A^^ 
20 (iV2 _ 1)4 



+ z^^ (-6208 + 49952 - 175780 N'' + 360896 iV^ - 469546 N" 
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+ 395168 N^^ - 205477 N^'^ + 56246 A^^^ - 6169 N^*^) 
1 

X 



24 {N^ - (iV2 - 1)^ 
+ z^^ (-204800 + 2013696 - 8684352 + 21407856 

- 33518156 + 34888595 iV^° - 24526859 N^^ + 11475108 N^^ 

- 3337802 N^^ + 545557 N^^ - 39943 7V^°) 

1 

X 



56 (7V2 _ 4)3 _ 1)7 

+ z^^ (-253184 + 2885888 - 14879968 N"^ + 46028112 A^^ 

- 95149021 + 138282732 - 144193790 N^^ + 108442897 A^^^ 

- 58120117 A^^^ + 21234008 A^^^ - 4926461 A^^° + 646061 N^^ - 36361 A^^^) 

+ 0{z''), (El) 



X 



8 - - if 



f- log Z - Z + i ^2 _ 1 ^3 _ I ^4 _ 1 5 _ 



1 + llAr2 



4 5" " 6(Ar2-i) 



7-4 + Ar2- ISAT^ ^ 11-81 Ar2 + Ar4 + 29 A^e 
+ 2z^ ^ hz^ 



+ 5z 



+ z 
+ z 
+ z 



9 {m - if 

10 19 - 277 A^^ + 127 A^^ - 319 A^^ 

10 {N^ - if 

11 -34 + 26 A^2 _ + 587 A^^ - 452 

11 {N^ - 1)^ 

12 132 + 1045 - 1813 A^^ + 2158 A^^ - 558 A^« - 115 N'^° - 69 A^^^ 



4 (Ar2 - 4) (A^2 _ 1)5 
+ z^^ (-1472 + 5408 A^^ - 24716 A^^ + 48156 A^^ - 97577 A^^ 

+ 118419 A^^° - 65911 A^^^ + 18318 A^^^ - 2614 A^^^) 
1 

X 



13 {N^ - 4f - if 
+ z^^ (5568 - 91344 + 393204 A^^ - 964687 A^^ + 1168570 A^^ 

- 845302 A^^° + 432359 A^^^ _ 133335 ^v^^ + 18275 A^^^ - 885 A^^^) 
1 



X 



7 (Ar2 - 4)^ {N^ - if 
+ z^^ (-25024 + 122896 A^^ - 380772 A^^ + 569163 A^^ - 1634619 A^^ 

+ 3326040 - 3347466 N^^ + 1923822 A^^^ - 642978 N^^ + 116764 N^^ 
- 7576 AT^O) 

1 



X 



15 (A^2 _ 4)3 ^^2 
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+ z^^ (43712 - 656912 N'^ + 3472772 - 10632331 

+ 18690504 - 20221284 N^^ + 15125616 A^^^ _ 3770834 A^^^ + 4069272 N^^ 
- 1342540 N^^ + 260364 N"^^ - 24083 A^^^) 

+ (E2) 



1 1 , l7^ 



16(Ar2-4)^(7V2_i)8j 

X = 1 + 3 ^ + 6 ^2 ^ 12 + 24 + 48 + 90 + 174 

o-28 + 29A^2 ^„ q 36-74 A^2^ 39 A^^ 
+ 12^^ + 18;s^ — -2 

m -203 + 625 - 619 + 232 A^^ 
(Ar2 - if 

„ 388 - 1592 A^2 + 2469 A^^ - 1725 A^^ + 469 A^^ 

^ 

736 - 3056 A^2 + 4302 A^^ - 3407 A^^ + 943 A^^ 



+ 6z 



12 



(A^2 _ 1) 



4 



+ 6z 



13 



1398 - 5850 Ar2 + 9318 A^^ - 6743 A^^ + 1922 A^^ 



{N^ - 1) 



4 



g -10520 + 47078 A^^ - 82728 A^^ + 70812 A^'' - 28891 A^^ + 3951 N^^ 

+ 6 z'^^ (79712 - 380208 A^^ + 730934 A^^ - 716200 

+ 372088 A^^ - 91009 N^^ + 8228 N^"^) 
I 



X 



(7V2 - Af {N^ - If 

+ 6 (600192 - 3641312 A^^ + 9388472 A^^ - 13344618 AT^ 

+ 11337734 - 5813592 Ar^° + 1719475 A^^^ - 267280 A^^^ + 16906 A^^^) 
1 

^ - Af - if 
+ 6 (4531200 - 33355776 A^^ + 107774336 A^^ 

- 200423616 A^^ + 236373604 A^^ - 183407604 Ar^° + 93969088 A^^^ 

- 31097639 A^^^ + 6341524 A^^^ - 722188 A^^^ + 35215 A^^°) 

1 

X 



(Ar2 - Af {N^ - if 
+ 6 z^^ (-8505600 + 71553024 A^^ - 268314976 A^^ 

+ 590267696 A^^ - 843497305 A^^ + 818280951 - 546354607 A^^ 
+ 249146769 A^^^ - 75461014 N^^ + 14404142 A^^^ - 1562138 N^^ 
+ 73433 Ar22) 



X 



(Ar2 - Af {N^ - if 
+ 6 z^^ (16021504 - 151615488 A^^ + 647496704 A^^ 
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- 1646428032 + 2769093400 - 3234820688 N^^ + 2679887383 N^^ 

- 1577573745 N^^ + 650694030 A^^^ - 182224525 N^^ + 32769509 N^'^ 

- 3399903 N'^^ + 154769 A^^^) 



1 

X 



+ 6 (-270452736 + 2603295744 - 11339998976 iV^ 

+ 29517395456 - 51057737736 A^* + 61725236828 iV^° - 53378649928 N^'^ 
+ 33218024089 iV^^ - 14765553520 A^^*^ + 4594385593 A^^^ - 966784918 A^^° 
+ 129788872 N^'^ - 9932849 N^^ + 325065 N^^) 

X ^- — 2 ^ + 0(z^^), (E3) 

..2 4 1 . 8 n + 1 „ 7-7-10Ar2 + 5Ar4 

G 3 3 iV2-l 2>{m-lf 

. -?,-7N^ + AN^ + N^ q -12-20Ar2 + 29A^^ -6Ar8 
+ 8z« — -3 + 8^^^ ,,,,, ,,4 



+ %z 



- ly 3 {N^ - ly 

10 12 - 6 A^2 + 15 A^^ + 8 A^6 _ n 

{m - if 

11 -82 - 34 Ar2 + 165 A^^ - 247 + 36 A^^ 



+ 8z'' 



+ 8z'' 



3{m - if 

-160 - 432 + 698 A^^ - 409 A^'^ - 2 iV^ - 18 A^^^ 
- 4) {m - if 

-1776 - 12392 A^^ + 3555 A^^ + 1850 A^^ - 1811 A^^ - 666 N^^ + 73 A^^^ 



3 (A^2 _ 4)2 ^^2 _ 1)4 

+ 8 z^"^ (3328 + 9920 A^^ + 23760 A^^ - 41452 A^^ + 12242 A^^ 

+ 5969 N^^ - 5077 A^^^ + 657 A^^'' - 5 A^^^) 
1 



X 



(A^2 _ 4)3 (^2 _ 1)5 

+ 8 z^^ (-114432 + 400128 A^^ - 1536160 A^^ + 1733488 A^*^ 

- 323935 - 602516 Ar^° + 449683 A^^^ - 124964 A^^^ + 5895 A^^^ 

+ 1675 N^^ - 75 N^°) 
1 

^ 3 - 4f (A^2 _ 1)6 
+ 8 (-91136 + 441856 A^^ - 1078144 A^'^ + 768256 A^'' 

+ 450812 A^* - 1054246 A^^° + 839396 A^^^ - 464465 iV^^ + 194981 A^^^ 

- 54586 N^^ + 9491 - 690 N'^^) 



X 



(Ar2 - Af (A^2 _ 1)7 
+ 8 z^^ (-425728 + 1766144 A^^ - 1985440 A^^ - 5221520 A^^ 
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+ 17290561 - 24613836 N^'^ + 24266712 N^^ - 17725908 N^'^ 
+ 9035120 N^'^ - 2989557 N^^ + 594069 A^^° - 63766 N^^ + 2193 A^^^) 
1 



X 



3 - 4)^ {m - if 

+ 8 z^^ (-1661184 + 3794944 + 7844000 

- 54955728 + 119244695 A^^ - 171667137 iV^° + 175851305 N^^ 

- 120836371 A^^^ + 53853632 A^^^ - 15114634 A^^^ + 2512157 N^^ 

- 229837 Ar22 + 9716 A^^^ - 526 A^^^) 

X I u 0(z^^) 

(Ar2 - 9) (A^2 _ 4)4 ^^2 _ 1)8 V ^ 

i„4-18A^2 + 32A^^-27A^6 + 18A^^ „ 
2 [m - 1)^ 



+ 6z'' 



3 - 16 Ar2 + 32 A^^ - 29 A^^ + 19 A^^ 



{N^ - 1)^ 

+ z^^ (49920 - 268032 A^^ + 595744 A^^ - 732080 A^'' 

+ 536515 A^^ - 239924 + 65162 A^^^ _ iQi^i N^^ + 829 A^^^) 
1 



X 



3 {N^ - Ay {N^ - ly 

+ 3 z^^ (92 - 890 + 3756 A^^ - 9052 A^^ + 13836 A^^ 
- 14023 A^^° + 9380 A^^^ _ 3590 tv^^ + 720 A^^^) 
1 

^ 2 (A^2 _ 1)8 

+ 3 z'^^ (36864 - 368640 A^^ + 1666560 A^^ - 4505344 A^^ 

+ 8111248 - 10267152 A^^° + 9373736 A^^^ - 6230112 A^^^ 

+ 3011024 A^^^ - 1020866 A^^^ + 227573 A^^° - 30248 N^^ + 1918 A^^^) 
1 



X 



{N^ - 4)^ {N^ - if 



+ z^^ (7215630336 - 96566722560 A^^ + 592839316224 A^^ 

- 2214635962752 A^^ + 5639017965384 A^^ - 10385973650808 A^^° 

+ 14322131449584 A^^2 - 15072994434276 A^^^ + 12212762091344 A^^^ 

- 7627059006688 A^^^ + 3656946091530 A^^° - 1335009872497 A^^^ 

+ 366664198242 A^^^ - 74601540139 A^^^ + 11013519148 AT^^ - 1144634371 A^ 
+ 79672354 A^^^ - 3373733 N^^ + 68030 A^^^) 
1 

X 



4 (A^2 _ 9)^ (;V2 - 4)^ (Ar2 - 1)^^ 
+ 3 z^^ (-275712 + 3201536 A^^ - 17022176 A^^ 

+ 54825776 A^^ - 119171333 A^^ + 184722213 A^^° - 210428660 A^^^ 



45 



+ 178609366 N^^ - 112928266 N^^ + 52388456 N^^ - 17133329 N 
+ 3699171 N^^ - 474540 N^^ + 29132 A^^^) 

X ^ q + 0(^28). 

{m - 4)^ {m -if ^ ^ 
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